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Summaky: A theory of creep buckling is presented in which the instantaneous 
elastic and plastic deformations following the application of a load, as well 
as the steady creep deformations, are considered in an approximate manner. 
Equations are given from which the critical time, that is, the time elapsing 
between load application and the collapse of the column, can be computed. 


1. Introduction 


Creep is the property of a solid to deform continuously with time under the 
action of a constant load. In the structural metals used in aircraft it is generally 
not observed at ordinary temperatures. It becomes noticeable when the temperature 
is increased, and it attains considerable importance at the temperatures to which 
aerodynamic heating raises structural elements exposed to the flow of air at Mach 
numbers equal to or greater than 3. 


Creep has an interesting consequence in elements subjected to compression, as 
can be easily explained with the aid of Fig. 1. The column shown there is 
assumed to be slightly curved in its natural state before the compressive load —P 
is applied to it. This is a reasonable assumption because a good mechanic can 
machine a bar to small tolerances, but never to zero tolerance. Similarly, an 
experienced testing machine operator can align the column with a great deal of 
accuracy, but never perfectly. The distance y between the line of action of the 
load —P and the centroid of the cross section of the column serves as a lever arm, 
and the product of this lever arm with the load is the bending moment to which 
the particular cross section is subjected. 


When the material is subject to creep, the bending moment increases the 
curvature of the axis of the column as time goes on. The increasing curvature 
causes increasing deflections y, and thus the bending moment acting on the section 
becomes larger. This vicious circle ends when the column snaps through and 
becomes useless as an element of the structure. 


*The work here reported was done at the Polytechnic Institute of Brooklyn under Contract 
No. AF33 (616)-116 with the Wright Air Development Center, Wright-Patterson Air Force 
Base, of the U.S. Air Force. The author is indebted to Dr. Sharad A. Patel of the Polytechnic 
Institute of Brooklyn for his help in the calculations. 
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Fig. 1. 
Column. 


Probably the most interesting fact about creep buckling is that it takes place 
under any compressive load, however small. One cannot talk therefore about a 
critical load in creep buckling. The significant quantity is the critical time, that is, 
the time that is needed for the column to collapse. The designer must know what 
the critical time is; he cannot provide a column which can stand up under the 
loads for an indefinitely long time; but he must make sure that buckling does not 
take place during the design lifetime of the structure. This design lifetime is often 
measured in minutes or even seconds in the case of missiles. 


Creep buckling was discussed in some detail by the author in an earlier 
publication”). There the deformations were assumed to be caused entirely by the 
secondary state of creep characterised by the straight line in Fig. 2. The 
instantaneous elastic and plastic deformations following the load application, as 
well as the transient creep phenomenon, were disregarded. 


In recent work done at the Polytechnic Institute of Brooklyn by Kempner, Patel, 
Erickson, and others'*’, the elastic deformations were found to be of considerable 
importance. Moreover, Odqvist showed in a recent note’? that integration of the 
differential equation of creep buckling is possible when the instantaneous plastic 
deformations and the transient creep are taken into account in an approximate 
manner. 
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Fig. 2. 
Representation of compressive creep strain-time curve when m=1 and n=3. 


The introduction of plastic deformations gives rise to considerable complica- 
tions in the analysis of creep buckling. As the plastic deformations are irreversible, 
the equations governing the deflections of the column must be modified when the 
compressive strain in the flange on the convex side of the column begins to decrease, 
and again when the compression in the flange changes into tension. It is the 
purpose of this paper to follow up all these changes and to present a theory of 
creep buckling when instantaneous elastic and plastic deformations take place 
and the material of the column is subject to transient and secondary creep. 


Research on creep buckling has been in progress in various institutions during 
the past few years. The many papers published on the subject are not quoted here 
since this paper is not based on them and since they are not needed for an under- 
standing of the analysis presented here. Reasonably complete lists of publications 
can be found in Refs. 1 and 8. A general discussion of the effect of creep on stress 
distribution and stability is given in Ref. 8. 


NOTATION 
a non-dimensional amplitude of first harmonic of w 


a, non-dimensional amplitude of first harmonic of initial deviations 
(amplitude divided by 4/2) after instantaneous loading when creep 
deformations begin 


a,, non-dimensional amplitude of first harmonic of initial deviations 
(amplitude divided by 4/2) in natural state of column before loading 


a, upper limit of non-dimensional deflection in integral 
a,,_ non-dimensional amplitude of second harmonic of w 


Q, non-dimensional deflection amplitude at buckling in first phase of 
creep deformations 
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non-dimensional deflection amplitude at buckling in second phase of 
creep deformations 


upper limit on a. OF Gor’ 

lower limit on @., OF dey’ 

value of a when strain reversal occurs in lower flange 
total cross-sectional area of two flanges of column 
strain quantity defined in equation (24a) 

strain quantity defined in equation (33) 

Fourier coefficient 

strain quantity defined in equation (245) 

Fourier coefficient 

Young’s modulus at elevated temperature 

distance between flanges of column 

constants 

constant in differential equation defined in equation (42) 
length of column 

exponent in creep law of equation (5) 

bending moment 

exponent in creep law of equation (5) 

axial force, positive when tensile 

radius of curvature 

time 

time required for buckling after tensile stresses appear in convex flange 
lower limit on critical time 

upper limit on critical time 

time required for tension to develop in convex flange 
simplified expression for fin, 

2y/h, non-dimensional total lateral deflection 


= 2y,/h, non-dimensional initial deviation after instantaneous loading 


when creep deformations begin 


= 2y,./h, non-dimensional initial deviation in natural state of column 


before instantaneous loading 

x’/L, non-dimensional axial co-ordinate 

axial co-ordinate 

total lateral deflection of centre line of column 


initial deviation of centre line of column from line of action of force 
after instantaneous loading when creep deformations begin 


initial deviation of centre line of column from line of action of force 
in natural state of column before instantaneous loading 
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€ strain, positive when tensile 


Strain in concave flange 
strain in convex flange 
= (h/L)’, Euler strain 
A material constant 
material constant 
radius of gyration of cross section of column 
o = P/A, average stress in column, positive when tensile (Ib. /in.*) 
, average stress in column when a=1 
7, stress in concave flange 
a, Stress in convex flange 


Tmax Maximum average compressive stress in column when loaded 
instantaneously 


2. Deformations 


For small deformations the curvature (1/r) of the centre line of a column can 
be expressed by the equation 


1 


where x’=xL is the longitudinal co-ordinate measured from one pin-supported end 
of the column and y=(h/2) w is the lateral displacement (see Fig. 1). The section 
is assumed to be an idealised I with the area concentrated in the flanges. The web 
has a negligibly small area, but is perfectly rigid in shear. The cross-sectional area 
of each flange is A/2. The rate of change of the curvature is 


( h ow 
where the dot indicates differentiation with respect to time 1. 


If the strain < is considered positive when tensile and suffixes “c” and “t” refer 
to the concave and convex flanges, respectively, the following geometric relationship 
holds, 


Differentiation with respect to time and substitution in equation (2) give 
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3. Stress-Strain Law 


The stress-strain law will be assumed as 


The stress o is positive when tensile. The first term in the right hand member /? 
represents the recoverable elastic deformations. The second term takes care of 
the instantaneous non-recoverable plastic deformations and, in an approximate 
manner, of transient creep. The quantities » and m are so chosen that the second 
term represents the intercept on the strain axis of the straight line of the secondary 
phase of creep, less the elastic deformation (see Fig. 2). The symbols m and n 
stand for integers and the value of k, depends on the conditions of loading : — > 


When « >0, *>0 (loading in tension): k,=1; 

When o>0, «<0 (unloading in tension): k,=0; 

When o <0, *+<0 (loading in compression): k,=1 provided that m 
is even; k,= - 1, provided that m is odd; 


When o <0, « >0 (unloading in compression): k,=0. 


Finally, the third term in the right hand member represents the non-recoverable 
creep deformations of the secondary phase. The values of k, are: — 


When 7 is even: k,=1 provided that o >0; k,=—1 provided that v <0; } 
When n is odd: k,=1. 


It is worth noting that the second term can also be written in the form 


m+ 
Ale ) 


The initial deviation of the centre line of the column from the straight line 
along which the load P is applied will be assumed in the form of a sine curve and 
will be designated as 


In a non-dimensional form 


W, = dy SiN TX. (6b) 


This assumption, and a similar one for the total non-dimensional deflections w at 
any time f, were shown to be satisfactory under certain conditions stated in Ref. 1. 
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When the column begins to creep, the compressive strain ¢. on the concave 
side of the column increases in consequence of the increasing curvature. For this 
flange of the column the deformation law is therefore 


when m and n are odd, as will be assumed for the rest of the calculations. At the 
same time, on the convex side of the column the increasing curvature diminishes 
the absolute value of the compressive strain <,._ Hence for this flange 


If the force P is taken positive when tensile and the bending moment M is 


(10) 


Substitution in equations (7) and (8), followed by a subtraction, yields 


If this expression is substituted in equation (4), the partial differential equation of 
the creep buckling process is obtained as 


4. Solution of the Differential Equation 


Following the procedure of Ref. 1, this equation can be reduced to a set of 
ordinary non-linear differential equations of the first order if the solution is 
assumed as 


and the various powers of w are expanded into Fourier series, Moreover, under 
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certain conditions a satisfactory approximate solution can be had if only the terms } 
multiplied by sin x are considered and terms containing the sine of a multiple 
angle are disregarded, as was shown in Ref. 1. 


This integration process will now be carried out under the assumption that 


In Ref. 7, Odquist made this assumption in his calculations. Moreover, incomplete 
experimental evidence available at the Polytechnic Institute of Brooklyn indicates 
that this choice of the exponents corresponds to the behaviour of 24S-T4 aluminium 
alloy bars at 600°F. 


The product ww =aa sin? 7x 


must be expanded into a Fourier sine series valid forO<x<1. The coefficient b, 
of the first term of the expansion is 


b,=2aa | sin*® txdx = aa... (15) 
0 
With n=3 the expression in square brackets in equation (12) becomes | 
—w)"—(1+w)"J= -2Bw+w%). . . (16) 
In the right hand member of this equation w* has to be expanded into a Fourier 
series. The coefficient c, of the first term of the expansion is 
) 
| sin' = (3) a... . (7) 
Substitution in equation (12) gives | 
This can also be written as \ 


a [26.42% - (2) (1+ Fa) ]=-(F) 


where ¢, is the Euler buckling strain 


(21) 
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In this non-linear ordinary differential equation, writing @=da/dt and 


separating the variables, 


% 
Integration yields 
ie tog [= -28 tan } (23) 
where A’=tpt+ - (2) ? . (24a) 


Equation (20) has a singular point at which 0a/ dt increases beyond all limits; 
naturally, the integration indicated in equation (22) should not be carried out beyond 
this point. The right hand member of equation (20) is always positive, while in 
the left hand member the expression in square brackets vanishes when 


If this value of a is less than or equal to unity, the lateral velocity of the column 
increases beyond all limits as the deflection approaches this value. If the value 
exceeds unity, equation (23) ceases to be valid, as shown in the next section. 


Naturally an indefinitely large lateral velocity means that the column buckles. 
Hence the time that elapses while the initial deviation amplitude a, increases to the 
amplitude a., can be designated as the critical time f.,. It can be computed by 
substituting the value of a., as obtained in equation (25) for the upper limit a in 
equation (23). 


5. Tension in Convex Flange 


The solution given in equation (23) is valid only as long as the stress in the 
flange on the convex side of the column is compressive for all values of x. It 
follows from equation (10) that the compression changes to tension in the middle of 
the column when a=1. The corresponding time is designated as the limit time fim. 
Replacement of a by 1 in equation (23) gives 


hin = E 2B tan (26) 
February 1956 9 
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When a, is small as compared to unity, this equation can be written in the 
simplified form 


Beginning with this limiting value of the time, the stress law of equation (8) 
must be replaced by 


in part of the column. In the remainder, equation (8) is still valid because the 
tension caused by the bending moment M is smaller than the compression caused by 
the force P. Solution of the problem under these conditions is difficult. Fortunately, 
the correct solution can be easily bracketed between upper and lower limits. The 
upper limit on the time necessary to reach a prescribed deflection is obtained if the 
plastic and the transient creep deformations are disregarded in the flange on the 
convex side of the column. Then equation (28) is not needed, and the solution 
given in equation (23) is valid at all times. 


This statement naturally is subject to the restriction that the upper limit of 
the integration in equation (22) must not be greater than a., in equation (25). If a. 
exceeds unity, substitution of its value into equation (23) yields an upper limit ¢.,., 
for the critical time. 


The lower limit of the time necessary to reach a prescribed displacement can be 
calculated on the assumption that equation (28) is valid for the entire length of the 
flange on the convex side of the column as soon as t> fim. In this case equation 
(11) must be replaced by 


t,-— w f=)" wy" w- (< ry’ (29) 


When m=1 and n=3, this equation becomes 


= - wiv-2(2 y (3w + w*). (30) 


Thus equation (20) is replaced by 
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The time ¢ necessary to increase the amplitude of the non-dimensional 
deflection from a, to a, is 


t= a(4+a) (44a) da (32) 
| 
” 
where A” (33) 
Integration gives 

1 log a,’ 4B tan 44 aa. (34) 


When the lower limit is unity and the upper limit a, the integral becomes 


The singularity of equation (31) corresponds to a displacement amplitude a.,’ 
given by 


tage =) 


Substitution of this value for a in equation (35) yields the time f,4 necessary for the 
lateral deflection amplitude a to increase from unity to a... Consequently, the 
lower limit ¢.., of the time required for buckling is 


6. Very Large Initial Deflections 


In some cases the deflection amplitude a, at the beginning of the creep process 
may be sufficiently large to cause tension on the convex side of the column. Then 
equation (29) characterises the deflections of the column during the entire creep 
buckling process provided that the simplifying assumption is made that the entire 
convex flange is in tension. The upper limit ¢.,, of the critical time is again 
obtainable through substitution of a,, from equation (25) into equation (23). The 
lower limit f,-; can be calculated by substituting a, for a, and a.,” for a, in 
equation (34). 


7. Deformations at the Moment of Load Application 


In the preceding analysis the quantity a, was the amplitude of the first sine 
component of the deviations of the centre line of the column from the straight line 
along which the load was acting at the moment when the steady creep deformations 
began. This quantity can be calculated with the aid of the laws of statics from 
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~ the amplitude a,, of the initial deviations which prevail before the load was applied. 
If the initial deformations before load application are designated by w,,, the 
increment in curvature is related to the strain in the flanges through the equation 


With the notation introduced in equation (19) the equation of equilibrium (equations 
(10) ) can be written as 


+w) 


(39) 
o,=o(1—w) 

During the loading process the absolute value of *=P/A is increased rapidly from 

zero to its final value; however, the loading is considered slow enough to render 

inertia effects negligible. As long as the compressive stress increases in both flanges, 

it is reasonable to assume that the instantaneous deformations are governed by 

the equation 


(40) 


Differentiation with respect to time of the expression in the right hand member 
of this equation gives the first two terms in equations (5) and (7) when m=1 and 
n=3. These terms represent the intercept on the <-axis of the straight line of the 
steady or secondary creep deformations (see Fig. 2). Substitutions and manipula- 
tions yield the differential equation 


where K*= (=) [ - + (<) |. (42) 


If the initial deflected shape is 


(41) 


Woo =Qoo SIN TX ; (43) 


the solution of the differential equation is 


w=a, sin x (44) 


with a,= ‘ . (45) 
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8. Strain Reversal 


It follows from equations (39) and (45) that on the convex side of the middle 
of the column (x=4) the compressive strain first increases, and later decreases. The 
reversal takes place when 


do, da, 
— —C 
dr 


do . . (46) 


Differentiation and algebraic manipulations yield 


[en 


Ayo = 


as the criterion of strain reversal. 


It can be seen from equations (45) and (47) that for any non-vanishing a,, 
strain reversal must take place before the denominator in the expression for a, 
vanishes. Consequently the column cannot buckle during the first phase of the 
loading process. 


The values the various quantities have at the moment of strain reversal will be 
indicated by a subscript “r.” From that time on the strain on the concave side is 
still governed by the pill 


but on the convex side the altered relationship begins to hold 


. (49) 


The range of validity of equation (49) spreads from the middle of the column as 
the load is increased. If the equation is assumed to be valid over the entire flange 
as soon as equation (47) is satisfied, the differential equation becomes 


h? ( ‘ 


Substitution of the earlier assumptions, expansion of (1+ w*) into a Fourier sine 
series, retention of only the first term of the series, and algebraic manipulations yield 
the following quadratic in the unknown a, :— 


This is the connection between a,,, a), and « beyond the strain reversal. 
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9. Buckling During Load Application 


If the curve has a horizontal tangent in a plot of the average stress o against 
the amplitude a, of the lateral deflection, the column buckles when the stress « 
corresponding to the horizontal tangent is exceeded. The slope @7/0a, is given by 
the implicit relationship 


In this equation the expression in square brackets is always negative. The right 
hand member of the equation is positive at the moment of strain reversal, as can be 
seen if the value of a, is substituted from equation (45) and subsequently that of 
a,, from equation (47): — 


(52a) 


The absolute value of the first term in the right hand member is certainly smaller 
than 4(o/,)?; hence the right hand member is greater than ¢, +(o0/E)-(c/,)’, 
which is known to be positive at strain reversal. Consequently, the slope 00 /0a, 
is negative and the compressive load in the column is still increasing at the moment 
when the strain reverses in the convex flange. 


When the absolute value of « is increased further, the right hand member of 
equation (52) decreases. This indicates that an increase in the deflections is 
accompanied by an increase in the compressive load (that is, by an increase in the 
absolute value of o, which is negative). The slope changes sign when 


This equation, in conjunction with equation (51), defines the values of a, and o at 
which the column can buckle during the second phase of the loading process. a, 
can be eliminated between this equation and equation (51). The following 
expression is obtained for the computation of the value of 7 at which buckling is 
possible during the second phase of the loading process: — 


The Aeronautical Quarterly 


10 


f 

re 

S 

| 
14 


ht 


of 


CREEP BUCKLING 


10. Tension in Flange During Load Application 


Equation (50) loses its validity when the stress on the convex side becomes 
positive. This happens when a,=1. The condition is (from equation (51)) 


When co is increased further, the strain on the convex side is 


The differential equation becomes 


alae ae ) = (1 —w,)?. (57) 


If the same procedure is followed as before, this reduces to the quadratic 


It is to be noted that in this equation the multiplier of (a /)° differs by 0-7 per cent. 
from that in equation (51) when a,=1. This inaccuracy is a consequence of having 
retained only the first terms in the Fourier expansions of the various functions. The 
slope of the curve of « against a, can be computed from the relationship 


The slope vanishes and the column buckles during the third phase of the loading 


process when 
3x 
Between equation (58) and this equation a, can be eliminated: — 


Equation (61) is the condition for buckling of the column during the third phase of 
the loading process. 
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11. Numerical Example 


A numerical example illustrates the application of the theory to practical 
problems. Results are given with slide rule accuracy. The material is characterised 
by the values 


m=1, n=s, E=7-4x 10° Ib. /in.’, 
u=0°8 x 10° Ib./in.2, A=87-5 x 10* Ib. in.~? hr.!/* 


Incomplete test results obtained at the Polytechnic Institute of Brooklyn indicate 
that these values can be used for extruded 24S-T4 bars at 600°F. 


The column has a slenderness ratio L/p=2L/h=50 and the amplitude of its 
half-sine shaped initial deviations from straightness is defined by the non- 
dimensional ratio a,,=2y,,/h=0:06. At the moment of load application, before 
creep can take place, the instantaneous deformations of the material caused by 
compressive stresses are characterised by equation (40). With the numerical values 
given, this relationship becomes 


= 13-S1o —0-782 x 10-‘e°. 
Equation (62) is plotted in Fig. 3. 


The behaviour of the column under rapidly applied loads, before creep 
deformations can develop, is shown in Fig. 4. The connection between the average 
compressive stress * and the non-dimensional deflection amplitude a, was calculated 
from equation (45) for the first phase of the loading. The stress o prevailing at 
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Fig. 4. 


Static deflections of column 
m=1 n=3 
L/p=S0 a, =0°06. 


strain reversal was obtained from equation (47) by an inverse method in which 
several stress values were assumed, the corresponding values of a,, were computed, 
and the value of o was interpolated at which strain reversal occurs in a column for 
which a,,=0-06. Upon substitution of «,, equation (45) gave the value of a, at 
strain reversal. The results were 


o,= — 17,910 Ib./in.’, a, =0°23)1. 


The connection between o and a, during the second phase of the loading is 
defined by equation (51). Solutions of this quadratic in a, were used in plotting the 
curve until the stress in the convex flange became tensile. This happens when 
a,=1. The corresponding stress is defined by equation (55) whose solution yielded 


= —22,320 Ib./in.’. 


Points on the third portion of the curve were computed from equation (58). The 
maximum average compressive stress during the loading process, which corresponds 
to static buckling, was obtained by a trial-and-error process from equation (61). Sub- 
stitution of this value in equation (60) yielded the corresponding value of the non- 
dimensional deflection amplitude. 
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The values obtained were 


Tmax= — 22,450 Ib./in.’, 


u = 11°67. 


teru=16°3 min. 
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History of mid-point deflections in creep buckling 


8 10 12 
t : MINUTES 
Fig. 5. 
n=3, 
oo = 0°06, a, =0-15, 
0°67. 
@ = 


14 


In the analysis of the creep buckling process it was assumed that the load was 
applied rapidly to the column until an average stress of 15,000 Ib./in.* was reached. 
Fig. 4 shows that at that time a,=0-15 in the absence of creep deformations. 
subsequently the load is left acting on the column and the material is allowed to 
creep in accordance with equation (5), the connection between time and deformation 
is furnished by equation (23). The points on the first part of the creep curve of 
Fig. 5 were computed from this equation. The first phase ends when the stress in 
the convex flange becomes tensile. This happens when a=1. Beyond that point the 
creep buckling process is characterised by equation (35). 


The column buckles when its lateral displacement velocity becomes infinite. 
From equation (25) the upper limit on a for infinite velocity is 


Substitution of this value in equation (23) gives the upper limit on the critical time. 


The Aeronautical Quarterly 


An 
pro 


Sub 


The 


N. J. HOFE 
| 
| 
| | 
| 
| | 
ie) 2 4 6 |_| 16 
pul 
16 
apy 
12 
de: 
int 
co 
mi 
va 
lo 
int 
is 
th 
Tl 
eq 
te! 
be 
st 
Be 
Fe 


CREEP BUCKLING 


The time necessary to reach a=1 is, from equation (26), 
tim= 11-95 min. 


An infinite lateral velocity is reached during the second phase of the buckling 
process when equation (36) is satisfied. Evaluation of the equation yields 


= 6°43. 


Substitution of this value in equation (35) gives the additional time 
ta=3°9 min. 
The lower limit on the critical time is, in agreement with equation (37), 
min. 


As the upper and lower limits are close, it is accurate enough for engineering 
purposes to state that the particular short column investigated buckles in about 
16 minutes after a load equal to two-thirds of the short-time buckling load is 
applied to it. By short-time buckling load is meant the critical load at the appro- 
priate temperature, not at room temperature. 


12. Conclusions 


A theory of creep buckling in which the instantaneous elastic and plastic 
deformations as well as the transient and secondary creep phenomena are taken 
into account has been presented. When the theory is used in the analysis of 
columns, the deflections of the column immediately following the load application 
must first be calculated. The increase in the deflection amplitude from the initial 
value a,, in the unloaded condition to the value a, prevailing immediately after the 
load is applied can be obtained from equation (45) when the compressive strain 
increases monotonically in both flanges during the load application. When this 
is not the case and when the strain on the convex side begins to decrease before 
the full value of the load is reached, a, must be computed from equation (51). 
The average stress + prevailing at the moment when the strain reverses is given by 
equation (47). The compression on the convex side of the column changes into 
tension when equation (55) is satisfied. Beyond this stage of the loading equation 
(58) must be used for the computation of the deflection amplitude a,. 


Naturally, the column may buckle during the load application if the stress o 
becomes sufficiently high. The condition of buckling is presented in equations (54) 
and (61). Equation (54) is valid when the critical condition is reached before the 
stress in the flange on the convex side becomes tensile, and equation (61) when 
the converse is true. 
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If the load acting on the column is not sufficiently high to cause buckling during 
the load application, the column will fail in creep buckling. Equation (25) gives 
the amplitude a., of the deflections at which the critical condition of infinite lateral 
velocity is reached provided that a.,< 1, that is the convex flange is still in com- 
pression. When the inequality is not satisfied, a,., still represents an approximate 
value characterising the critical condition. Substitution of this value for a in 
equation (23) yields the critical time necessary for buckling. The time is rigorously 
correct when a,, <1 and is an upper limit when the inequality is not satisfied. 


The time fin at which tensile stresses appear in the flange on the convex side 
during the creep process can be calculated from equations (26) and (27). If the 
critical deflection amplitude is reached beyond this time, it must be calculated from 
equation (36). Substitution of this value in equation (35) and consideration of 
equation (37) yield a lower limit for the critical time under these conditions. 


Finally, the deflections immediately following the load application may be 
sufficiently large to cause tension on the convex side of the column. Under these 
conditions the upper limit ¢.,, of the critical time is again obtainable through 
substitution of a,, from equation (25) into equation (23). The lower limit t...,; can be 
calculated by substituting a, for a, and a.,’ for a, in equation (34). 


REFERENCES 


1. Horr, N. J. Buckling and Stability, Forty-First Wilbur Wright Memorial Lecture, Section 
5.2, Columns Subject to Creep. Journal of the Royal Aeronautical Society, Vol. 58, No. 
517, p. 30, January 1954. 


te 


KEMPNER, JOSEPH. Creep Bending and Buckling of Linear Visco-Elastic Columns. 
N.A.C.A. T.N. 3136, January 1954; also P.I.B.A.L. Report No. 195, March 1952. 


3. KEMPNER, JosEPH. Creep Bending and Buckling of Non-Linear Visco-Elastic Columns. 
N.A.C.A. T.N. 3137, January 1954; also P.I.B.A.L. Report No. 200, May 1952. 


4. KEMPNER, JOSEPH and PATEL, SHARAD A. Creep Buckling of Columns. N.A.C.A. T.N. 
3138, January 1954; also P.I.B.A.L. Report No. 205, November 1952. 


5. Ness, NATHAN. Time Dependent Buckling of a Uniformly Heated Column. N.A.C.A. 
T.N. 3139, January 1954; also P.I.B.A.L. Report No. 192, October 1951. 


6. PATEL, SHARAD A.; BLOOM, MARTIN; ERICKSON, BURTON; CHWICK, ALEXANDER and Horr, 
N. J. Development of Equipment and of Experimental Techniques for Column Creep 
Tests. N.A.C.A. T.N. 3493, September 1955; also P.I.B.A.L. Report No. 239, December 
1953. 


7. Opgvist, Forke K. G. Influence of Primary Creep on Column Buckling. Journal of 
Applied Mechanics, Vol. 21, No. 3, p. 295, September 1954. 


8. Horr, H. J. Rapid Creep in Structures. Journal of the Aeronautical Sciences, Vol. 22. 
No. 10, p. 661, October 1955. 


20 The Aeronautical Quarter!’ 


| 

1. 

are 
| 
it 
co 
inc 
| 
| 
the 

be 
fo 
| 
W 
co 
of 
at 
Re 
Feb 


ns. 


Supersonic Bangs 


Part [ 
P. SAMBASIVA RAO 


(Department of Mathematics, The University, Manchester) 


SumMarRyY: In this paper a non-linear theory of the supersonic bangs from 
an accelerating body is worked out, using an extension of Whitham’s” 
method. The strength of the bow shock is determined, and it is shown that 
the effect of acceleration on the pressure rise is appreciable at large distances 
from the nose for Mach numbers near unity. The results are compared with 
those of linear theory; the estimate for the pressure rise is found to be about 
a half of the linear value. The results also include an expression for the 
curvature of an attached shock at the nose, and a formula for the decay of a 
shock after it is detached from a decelerating body. 


1. Introduction 


When an aircraft is flying at supersonic speed, sharp, loud and explosive sounds 
are heard by observers on the ground. These sounds are called “ supersonic bangs.” 
They are heard when the observer receives shock waves produced by the aircraft as 
it pushes through the air. There has been some discussion (Yates®, Warren’) 
as to whether the accumulated noise from the engine and other parts of the aircraft 
could contribute to the production of bangs, but this so-called “acoustic explana- 
tion” was rejected since it could not give the correct magnitude for the pressure 
increase, and, as Warren pointed out, bangs are heard from rocket-boosted missiles 
when the motor is no longer burning. 


For a body moving with constant supersonic velocity, the strength and position 
of shocks are predicted by the theory of Whitham"). In the present paper, that 
theory is extended to the case of accelerated motion. The theory is developed for 
a uniform atmosphere, but the method may be extended to other cases. 


The general features of the shock wave patterns for steady supersonic flow have 
been known for a considerable time. But it is only recently, since bangs were 
actually heard, that much consideration has been given to the shock wave patterns 
for accelerated motion or for motion which is not in a straight line. The earliest 
criterion to determine when a bang will be heard seems to have been given by 
Warren’, He has shown that a point receives a bang if at any previous time the 
component of the velocity of the aircraft towards the point was equal to the speed 
of sound. Gold“) has given the following simple form of the argument. In general, 
at an instant ft, a point receives a wavelet generated by the source (aircraft) at some 
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single instant 7; but under certain conditions it receives at the same time wavelets 
generated at successive instants. In the latter case, the wavelets combine to give 
a sudden increase of pressure at the point. Then a bang is heard at the point. 
Thus, if R is the distance of the source from a point at time 7, R+5R the distance 
at time 7+6r and a, is the speed of sound in the air, a bang is heard at the point 
at time ¢, if 


that is, =-a. . ‘ (1) 


This is the condition that the component of the velocity of the aircraft towards the 
observer is equal to the speed of sound. Gold has pointed out that an observer 
will most likely hear two bangs, one from when the source is accelerating and the 
other from when it is decelerating through the speed of sound, and the latter one 
will be the louder because it is generated when the aircraft is nearer. It will also 
arrive at the observer sooner. Of course, there will be more bangs if the component 
of the velocity of the aircraft towards the observer becomes sonic more often. 
According to this approximate method of determining the bangs, the points which 
receive a bang at time ¢ lie on the surface given by 


=7+(R/a,) 


and this is the envelope of the sound waves sent out by the body. Such envelopes 
correspond, in the simple acoustic theory, to the real shock waves; in this paper 
they will be referred to as “ wave fronts.” 


Lilley, Westley, Yates and Busing’ have determined the pattern of wave 
fronts due to a source accelerating and decelerating through the speed of sound 
along a straight line and moving with constant speed along certain curved paths. 
They have pointed out that a novel feature of the wave fronts in accelerated flow is 
the formation of cusps on the wave fronts and they have determined their positions. 
A typical case is shown in Fig. 1. They have checked the results experimentally 
with the wave patterns produced on the surface of shallow water when a body moves 
through it. (This makes use of the “hydraulic analogy,” i.e. the correspondence 
between the shallow water theory and the theory of compressible fluid flow). For a 
finite body, there is also a shock generated near the tail of the body (a finite body 
is like a doublet) and the actual flow pattern is as shown in Fig. 2. 


Reference 5 gives an approximate description of the geometry of shocks but 
gives no method for predicting their strengths. The pressure increase to be expected 
has been considered on a linear theory by Warren”, who obtains a formula for 
the maximum pressure increase at any point. He has checked his estimate for 
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2): 


(1) 
Fig. 1. Fig. 2. 
(1) Front portion of the wave front. (1) Body. (2) Front shock. 
(2) Cusp. (3) Rear portion. (3) Tail shock. (4) Wake. 


pressure rise in steady supersonic flow against the results of the more exact theory 
given by Whitham“? and has found that it is considerably greater than the correct 
value (in one case, by a factor of 13). Later, in an addendum to his paper, Warren“? 
discusses how the linear theory fits in with the true situation in which the shock 
waves caused by the motion of the aircraft are really responsible for the bangs. 


The inadequacies of the linear theory have been pointed out by Whitham”). 
In the linear theory, the disturbances are propagated with constant speed equal to 
that of sound in the undisturbed fluid. The variation of the local speed of sound 
and the convection of sound with the fluid are not taken into account. These 
errors accumulate as the disturbances propagate outwards and, to get a valid 
description at large distances, the appropriate speed of propagation (equal to the 
local speed of sound plus the local fluid velocity) must be taken. Again, the linear 
theory cannot give a true description of the shocks since they have essential non- 
linear properties. They are represented inadequately by the wave fronts of linear 
theory. For example, at these wave fronts the pressure is in fact continuous, 
increasing behind the wave front as the square root of the distance from the 
wave front. To remedy this, it is supposed that the linear theory gives a correct 
first approximation to the variation with distance of the values of pressure and 
velocity carried by a wavelet, but that its position requires correction. This gives a 
valid description of the flow at large distances, and shocks become a necessity to 
separate regions where otherwise wavelets would overlap. To a first approximation 
the speed of the shock at a point is the mean of the speeds of the two wavelets which 
meet at that point. From this condition the shocks may be determined. 


In this paper the linear theory of supersonic bangs is modified in accordance 
with the method just described. The effect of acceleration on the bow shock is 
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studied and the strength of the bow shock is obtained. The results of this theory, 
for constant velocity, reduce to those of Whitham". A comparison with the steady 
flow problem has shown that the influence of acceleration on the strength of the 
shock becomes appreciable at large distances from the nose for Mach numbers 
near unity, when the rise of pressure is found to be about a half of that given by 
the linear theory of Warren’. When the shock is attached to the body, the strength 
of the shock at the nose and the angle of the shock at the nose are independent of 
the acceleration of the body and therefore are exactly the same as the results of the 


steady flow problem. But the present theory shows how the shocks vary with 
acceleration, away from the nose. 


As already explained, the theory is obtained by a suitable modification of the 
linear solution to the problem. Now, although the full linear solution can be 
obtained, it is too complicated for the modification to be applied as it stands. 
Therefore, the approximate form which it takes near the wave fronts (which are 
ultimately replaced by shocks) is considered in detail. But the behaviour near wave 
fronts is precisely the subject of “ geometrical acoustics,” and valuable information 
can be obtained directly without going to the full linear theory first. Part of the 
information given by geometrical acoustics is the geometry of the wave fronts 
(discussed by Lilley and others‘). But, in addition, and this does not seem to be 
generally known, the approximate variation of flow quantities immediately behind 
the wave front can be determined. These results are conveniently found by assum- 
ing that the pressure increase (for example) can be expanded in the form 


(x, y, 2) F (t— (x, y, Z)) +x, (xy, +...... (3) 


where t— (x, y, z)=0 is the wave front and 


F,(¢)= F (A) dA, F,(c)= F, (A) dA, (4) 


At this stage, the function F (~) could be taken as a power of o (for flow past a 
pointed axi-symmetric body the appropriate power is c+); but it is taken in this 
more general form because it will be seen that in certain circumstances this expansion 
is also valid when (t—w) is not small. An approximate solution x (x, y, z)F (t—©) 
might be expected on physical grounds, since it represents wavelets given by 
t—w (x, y, z)=constant and with a varying amplitude given by x(x, y, z). When 
this expansion is substituted in the wave equation and the coefficients of F”, F’. 
F ...on both sides of the equation are compared, equations for , x, X;,.. . are 
obtained. The equation for w is the “eikonal equation” of geometrical optics and 
governs the geometry of the wave fronts. The equation for x then determines the 
variation of p in the neighbourhood of the wave front. The details are given in 
Section 2. 


When the solution can be approximated in the form x(x, y, z)F(t— ©), the 
correction can be applied quite straightforwardly. The approximate form is 
replaced by x (x, y, z)F (n), where (t, x, y, z) is determined such that the equation 
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»= constant represents the propagation of a wavelet more accurately than (t—)=0 
of the linear theory, i.e. an improved value for the propagation speed is introduced 
in place of constant value a,. 


Now a crucial point must be made. In the linear theory, geometrical acoustics 
(in the general form outlined) determines the flow quantities near the wave front. 
But this is not sufficient for the purpose of determining the shocks in the modified 
theory. For, in the non-linear theory, wavelets are continually running in to the 
shock as it propagates (this is one of the essential properties which is lost in the 
linear theory). Thus, in the improved theory, the value of » for the wavelets just 
behind the shock increases as the shock propagates out, since the earlier ones have 
been fed into the shock. Hence, to describe the shock at large distances from the 
source, the solution when » is not small is required, and this means that the linear 
solution is also required when (t—w) is not small. This is obtained from the full 
linear theory which is discussed in Section 3. It is found that the flow quantities 
may be approximated by expansions of the type of equation (3). The amplitude 
factor x is identical with the one obtained by geometrical acoustics, but now the 
region of validity of the expansion is extended and the function F is determined. 
This function is found to be exactly the same as the one which appears in the 
Whitham theory for the steady flow problem. (In that case x is inversely propor- 
tional to the square root of the distance from the flight path). 


The linear theory is presented in Section 3 for any axi-symmetric slender body 
moving along a straight line. In his paper, Warren considers a special shape for 
the body, but the general case is not more difficult. The details of the improvement 
are set out in Section 4. The main result is the expression for the pressure 
increase at a point as the shock passes through it. The results reduce to those for 
the steady problem when the acceleration is set equal to zero. In Section 5 the 
results of the theory are discussed in detail. First the behaviour of the shock near 
the nose is considered. As already noted, the angle and strength of the shock at 
the nose are the same as in the case of constant supersonic velocity, but the formula 
for the curvature of the attached shock at the nose is new; it is compared with 
the expression found by Cabannes‘’. Numerical results are then given which 
exhibit the effect of both acceleration and retardation on the shock strength at vary- 
ing distances from the aircraft. Finally, the theory is used to study the decay of a 
shock at large distances after it is detached from a body decelerating through sonic 
speed. 


In this paper, only the bow shock is studied. It is expected that, as in the 
steady flow problem", the shocks starting near the tail of the body (see Fig. 2) 
will be of approximately the same strength as the bow shock. When the acceleration 
is positive, as already pointed out, the wave front is cusped and there is a rear 
portion, as shown in Fig. 1. Only the front portion away from the cusp is 
considered in this paper. But the rear portion is of much smaller strength, since it 
is produced earlier, when the body is moving with smaller speed; hence it is less 
important. An important addition to the theory is to extend its scope to 
include curved flight paths. There appears to be no essential difficulty to this 
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and some results in this direction will be published in Part II. The 
main difference will be that the amplitude function x will depend on an angular 
co-ordinate about the direction of motion, and the shocks will be stronger on the 
inside of the curve than on the outside. 


NOTATION 


u,v 


Ss 
f (7), 8 (7), (7) 
X (7) 


R 


w (x, y, Z) 

A,o 

a 

Ay» Po» Po 


Cartesian co-ordinates of a point at time ¢ 
distances of a point along and perpendicular to the flight path 


components of velocity along and perpendicular to the flight 
path 


distance of any point along the ray 

Cartesian co-ordinates of the source at time 7 

distance of the source along the flight path 

time when the source crosses the ray through the point (x, r) 
distance between the source and the point 

distance of an element of the body from its nose 

maximum value of € 

values corresponding to £=&, 

accurate wavelet variable 


function related to the wave front 


functions occurring in the expansion of a physical quantity near 
the wave front 


variables 

speed of sound in the fluid 

pressure and density in the fluid 

corresponding quantities in the undisturbed fluid 
speed of the source 

Mach number (U/a,) 

(M? — 

Mach angle, sin-? (a,/U) 


component of acceleration of the source towards the point 
(x, y, Z) 


velocity potential 

Laplace operator 

area of cross section of the body 

direction cosines of the outward-drawn normal to the ray tube 
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m Strength of a source 
y ratio of specific heats of a fluid 
K 
b(s) defined by equation (51) 
€ semi-vertical angle of the nose of the body 
6 angle of inclination of the shock with the flight path 
k curvature of the shock at the nose 
k, curvature of the body at the nose 
5 a constant 
root of the equation F (y)=0 
d__ distance by which the detached shock is ahead of the body 


n, the characteristic on which the shock forms first 


2. Geometrical Acoustics 


Consider a point source moving with supersonic velocity in an otherwise 
undisturbed fluid, and suppose its position at an instant 7 is given by the co-ordinates 
[f(7), g(7), h(7)]. The disturbance caused by the source at time = will spread out, 
according to the linear theory, with speed a,, the speed of sound in the undisturbed 
fluid. At time ¢ this disturbance will reach the surface of the sphere 


If the spheres of disturbance generated by the source at different times have an 
envelope at time f, it is given by (5) and 


This is the wave front, the linear approximation to the shock. From relations (5) 
and (6) it is seen that 
— 3, =a, 
which is the condition that at time = the component of the velocity of the source 
towards the point (x, y, z) is equal to the speed of sound. This is the condition that 
the point (x, y, z) receives a bang and, from (5), the bang is heard at the point at 
time (R/a,) after the time r. Thus the wave front moves out at speed a, along the 
line drawn in the direction in which the component of the velocity of the source 
at the instant 7 is a,. Such lines are called “rays.” (Of course, there will be a 
whole cone of such rays from any point of the flight path; but, in view of the axial 
symmetry, discussion can be confined to one meridian plane). 
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A ray makes an angle (4 — ») with the flight path, where » is the instantaneous 
Mach angle sin-'(a,/U), and U is the instantaneous speed. The — of the 
wave front at time ¢ can be written in the form 


Then since, from (7), a,~'(0R/ér)+1=0, the direction ratios w,, ©,, ©, of the 
normal to the wave front are proportional to (0R/dx), (@R/dy), Hence 
the normal is parallel to the ray through (x, y, z). Thus the rays are the orthogonal 
trajectories of the wave fronts. 


If the wave front is cusped, the cusp satisfies relations (5) and (6) and 


which is obtained by differentiating (6) with respect to 7. From this the expression 
for the distance of the cusp from the point (f, g, h) is found to be 


where U =(f? + 9? +h’)! is the speed of the source and A is the component of the 
acceleration of the source towards the point (x, y, z). Since U >a, it is seen from 
(10) that the wave front will be cusped only if A is positive. The wave front 
for uniform positive acceleration is sketched in Fig. 1. 


Now the behaviour of the flow near the wave front is considered. It will be 
convenient to work with the velocity potential ¢; the pressure increase, which is 
given by p— p,= —/po%:, and the velocity components may then be deduced. In the 
linear theory ¢ satisfies the wave equation 


As already mentioned in Section 1, the solution is now expanded in the form 


where again f, [ f(A)dA, f.(7)= | 
0 0 

and ©, x, X; X2,--.. are functions of (x, y, z) only. Substituting this expansion for 


¢ in the wave equation and equating the coefficients of f”, f’, f, we get 


a, (w,? + 0,?+0,’)=1, 
2 (Mr + Oy XiytO2Xid +X, V? x, ° (15) 


Equation (13) is the “eikonal” equation of optics and it is satisfied by the wave 
front given by equation (8). 
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Equation (14) can be written in the form 


)+ dy (a, )+ oz (a, )=0, (16) 


from which the variation of x along a ray is deduced. If equation (16) is integrated 
over the domain enclosed by a ray tube and bounded by the wave front at two 
instants, then by the divergence theorem 


s 


where S is the surface enclosing the domain and (A, y, v) are the direction cosines of 
the outward-drawn normal at any point of the surface. Since ,, o,, , are the 
direction ratios of the ray through (x, y, z) there is no contribution to the integral 
in (17) from the sides of the ray tube. On the ends of the tube, where (A, », v)= 
+ (4,2, A, My, Az), Ay (AW, + wo, +vo,)= +1. Therefore, if AS is the cross-sectional 
area of the tube, (17) reduces to the difference between the values of x?AS on the 
two ends. Thus, along a ray, 


It may be noted immediately that as a cusp is approached the rays converge and 
AS tends to zero; then, from (18), x tends to infinity. 


Fig. 3. 


The amplitude x will now be determined, for a straight flight path, in terms of 
the velocity and acceleration of the body and the distance along the ray. If A and B 
in Fig. 3 are two positions close together on the flight path, BB, and AC, are the 
rays from B and A, and BC and B,C, are sections of the wave front, then, since 
the angle CAB is (4—-p), 


where R=AC,. Equation (18) then gives 


27x?RB,C, cos »=constant. ; j (20) 
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But if the increments in » and 7 from A to B are 5 and 87 respectively 


and AB=U5-. 
dp 
Therefore xo RU COs? (144 U cosa dr 


It may be noticed from (22) that x becomes infinite for R=A~-!(U?—a,”), where 
A denotes the component of the acceleration of the source along the ray; this agrees 
with the position of the cusp given by (10). 


Since p — p= — Po, it is found from relation (11) that the pressure increase is 
proportional to x near the wave front; hence (22) shows how the pressure increase 
in the neighbourhood of the wave front varies with R. Since F (7) is proportional to 


ot for small o, the pressure is continuous across the wave front, but dp/dt 
is infinite. 


As pointed out in Section 1, the solution given by geometrical acoustics for 
small o is not sufficient for the treatment of shocks in the “improved theory.” The 
solution for large R also is needed and this is obtained from a full linear theory. 
In particular the function F (~) is determined and, as may be expected, it turns out 
to be the same as the corresponding function in the steady flow problem. 


3. Linear Theory 


In this section, the linear theory is considered for the case of a slender, axi- 
symmetric body which moves in a straight line, through the fluid. The ensuing 
flow, assumed to be irrotational, can be described by a velocity potential ¢ (x, r), 
where x is measured along the flight path in the negative direction and r is the 
distance from the flight path. The perturbation velocity components u and v, along 
and perpendicular to the x-axis respectively are 


If the second and higher powers of (u/a,), (v/a,) and their derivatives are neglected, 
¢ satisfies 


and from Bernoulli’s equation the pressure increase is given by 


where p, and p, are pressure and density respectively in the undisturbed fluid. 
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Fig. 4. 


As the body moves, it pushes aside the fluid, and its effect at any time may be 
represented by a source distribution of appropriate strength along that part of the 
flight path occupied by the body at that time. In fact, as the body moves through 
a point (X, 0), the rate at which the fluid is pushed aside is US’ (£) per unit length 
of the flight path, where U is the speed and S(€) is the area of cross section of the 
body at a distance € from the nose. To obtain the potential at a point (x, r) at 
time f¢, it is necessary to sum over all such points (X, 0), that contribute to the 
disturbance at (x, r) at time ¢. The time 7 at which a disturbance must leave (X, 0) 
is given by 


If the position of the nose of the body at time 7 is given by —f(z), then the speed 
U=f (r), and the element of the body at a distance € from the nose is at a point 
(X, 0) given by 


The contribution of a source of strength m at (X, 0) to the potential at (x, r) is 


ial 
Thus, the potential is +r} j 


where the integration is over all points (X, 0) for which the source strength is non- 
zero at the time 7 given by (26). 
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To deal with the integration the geometrical representation in the X, 7-plane 
will be used. Equation (26) represents a hyperbola with its centre at (x, t) and 
semi-axes r and (r/a,). The integration is over the part of the hyperbola which 
gives values of € (in equation (27)) corresponding to the points on the body. The 
body is effectively infinite if the wake boundary is included and ¢ varies from 0 
to 00. The integration is over the part of the hyperbola for which X >-f (7). Of 
course, the main contribution comes from a range of € which is of the order of the 
length of the body. The case in which the body is moving throughout with subsonic 
speed is shown in Fig. 4. The curve CD is the trace of the nose and A,B,, A,B, 
are typical hyperbolae (26). For a given point (x, r), as ¢t increases the hyperbola 
moves upwards; thus A,B, corresponds to a greater value of ¢ than A,B,. The 
region of integration is the part of the hyperbola to the right of the path of the nose. 
(The upper branch of the hyperbola refers to disturbances for which 7 >t and need 
not be considered). For any value of t, the hyperbola intersects the flight path pro- 
vided that the body has been flying for a sufficiently long time. Thus, in purely 
subsonic flight, a point (x, r) will receive disturbances at all times. 


Figure 5 represents the case of an accelerating body having partly subsonic and 
partly supersonic flight. If ¢ is chosen sufficiently small, the hyperbola intersects 
only the subsonic part of the flight path, but not the supersonic part; this is so for 
A,B,. As t increases, a value is reached for which the hyperbola A,B, just touches 
the supersonic part of the flight path of the nose while intersecting the subsonic 
part. This corresponds to the arrival of the front portion of the wave front 
(see Fig. 1) at the point (x, r). As ¢ increases further the hyperbola penetrates 
into the supersonic part of the flight path as shown by A,B,. The integration in 
(29) is then over two regions; one contribution is from the hyperbola between A, 
and B, in the supersonic part of the flight path, and the other from the region 
beyond D in the subsonic part. Eventually there will be a second singularity in 
the expression for ¢ when the hyperbola moves up to touch the flight path between 
B, and D. This corresponds to the arrival of the second wave front, i.e. the rear 
portion shown in Fig. 1. When the two points of contact with the flight path 
coincide to give a triple point of contact, the two wave fronts coalesce to give the 
cusp. For a given ray, the wave front will always correspond to the hyperbola 
touching the flight path at the same point (X, r). As the wave front moves out along 
the ray the curvature of the corresponding hyperbola decreases. Initially it is greater 
than the curvature of the flight path and it becomes equal to that of the path when 
the wave front has its cusp on the ray; then it becomes less than that of the path as 
the corresponding point (x, r) moves from the front to the rear portion of the 
wave front. 


For any motion of the body the bow wave front will correspond, as in Fig. 5, 
to the hyperbola just touching the flight path, and for points close behind the wave 
front the hyperbola will have a finite section like A,B,. The solution ¢ is now 
considered for such points. There will, of course, be contributions from other 
sections of the hyperbola inside the flight path, such as the subsonic part beyond D 
in Fig. 5. But, provided that the neighbourhoods of cusps (where contributions 
merge) are excluded, these will, in general, be of smaller order. For they will have 
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Fig. 5. 


been generated earlier and travelled much larger distances; for example, in the case of 
Fig. 5, they are generated during the subsonic motion. In any case, these other 
contributions change relatively slowly and the sudden increase, which is crucial for 
describing the shock, is the contribution from A,B,. 


To evaluate the integral in (29) it is convenient to change the variable of inte- 
gration to €. In Fig. 5 for the hyperbola A,B,, there is a maximum value of € in 
the supersonic part A,B, of the hyperbola; let it be denoted by €,. This €, is the 
length of the body contributing to the supersonic part of the disturbance. If 
(X,, 7,) is the point of the hyperbola corresponding to £=€,, and R,?=(x—X,)? +r’, 
it will now be supposed that (€,/R,)«1. This approximation is sufficient for the 
determination of shocks and is similar to the one made in the problem of steady 
flow; it is a precise statement of the condition discussed in Section 1 for the deter- 
mination of the shocks at large distances. Since in any practical case the acceleration 
cannot be very large, it is reasonable to assume that the change in the speed in the 
time the body takes to fly its own length is small compared to its speed. That is, 
(IU’/U*) « 1, where | is the length of the body. Expanding € in powers of (X — X,), 
near €=€,, we have 

2 


neglecting higher order terms. This is equivalent, in the X, 7-plane, to replacing 
the hyperbola in the neighbourhood of =, by its circle of curvature at £=€,, which 
is permitted because of the approximation (£,/R,)«1 and the assumption about 
the acceleration. From (30), for € near &,, 
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In the integral (29), consistent with (€,/R,)«1 and (JU’/U?)«1, R and U(r) can 
be replaced by R, and U (r,) respectively. Therefore the approximate form for ¢ is 


S’ (6) df 
TR 2 x-1, 


From (26) and (27), for any &, 


dé U x-X 

For £=£,, (d€/dX) vanishes and therefore from (33) 

2¢ a, 
RY 
is 


The co-ordinates of the point (x, r) will now be chosen in a convenient way, based 
on rays. Let the nose of the body cross the ray through the point (x, r) at time 7, 
and let s be the distance of (x, r) along the ray from the nose. Then, taking (7,, s) 
as the new co-ordinates, the point (x, r) is 


on 


r= 


As the new co-ordinate in place of the time it is convenient to take ¢,, which is 
essentially a measure of the time which has elapsed at any point after the passage of 
the wave front. After the wave front has gone beyond the point (7,, s), the 
disturbance generated by the element of the body distant €, from the nose at time r, 
reaches the point at time t. Hence £, is approximately equal to the distance travelled 
by the body in the time (7,—7,). Since t—7,=s/a,, 
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Now, in time (7,—7,), the body moves a distance which is of the order of its own 
length and the changes in the values of U, U’ and & in this time are small and have 
already been neglected. In (35), 7, may be replaced by 7, and then it becomes 


—U?(r,) 


(38) 
T 2 1 al 


This gives the value of ¢ at a distance s along the ray, labelled 7,, at a time 
(£,/U (r,)) after the wave front passed through the point. 


Differentiating (38) and neglecting terms of order £,/s we have, in view of (37), 


(€,) [@,MB?s — M’s?]-3, 39 
— (€,) [a,MB?s — M’s?]-3, ) 
S” (6) dé 
where F(¢,)= 
U (7) M, U (7) M’, 
a, a 


¢, is the component of the particle velocity along the ray and it can be shown that 
the component tangential to the wave front is of the smaller order £,¢,/s. Hence 
the particle velocity is approximately normal to the wave front. The tangential 
component must, of course, be zero at the wave front in order to satisfy the law of 
conservation of tangential momentum across the wave front. From (39) it is 
seen that 


and it will be noted later that this again is a requirement of the conservation laws. 
It should be noticed that (39) is in the form suggested by geometrical acoustics 
in Section 2 and the amplitude is exactly the same. For M’=0, equations (39) 
reduce to those of steady flow". 


It may be mentioned again that the formulae apply only in the region near the 
front portion of the bow wave front. If the acceleration is positive, as already 
pointed out, the bow wave front is cusped and it will have a rear portion. The 
approximations made cease to be valid as the cusp is approached and the neighbour- 
hood of the cusp cannot be dealt with by this method. 
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4. Improvement of the Linear Theory 


The linear theory, given in Section 3, will now be improved, using the 
technique of Whitham"). More precisely, an extended version of that theory is 
required. The steady flow problem could be described by two independent 
variables, but in the unsteady flow, studied here, there are three independent 
variables (t, 7, s). (The subscript of 7, is dropped hereafter). Whitham has extended 
his method to include such cases and a detailed account of the technique with 
applications to other problems is to be published shortly. Briefly, it turns out that 
the propagation of wavelets along each ray may be considered separately. This is 
because only the derivatives ¢, and ¢, are important and 7; plays the role of a 
mere parameter. Along a ray, there are only two effective variables t and s, and 
the theory is rather like a one-dimensional unsteady flow outwards along each ray, 
with amplitude given by x. 


In the linear theory all disturbances propagate with constant speed a,, while in 
reality they propagate with the local speed of sound relative to the fluid. The 
technique accepts the variation in amplitude predicted by the linear theory but 
corrects the speed of propagation. When the correct speed is introduced, as 
explained in Section 1, the shocks can be included and their position and strength 
can be determined. 


As the wave front, which is the counterpart of the shock wave in the linear 
theory, advances into the fluid at rest, it is the normal component of the velocity 
that changes, while the tangential component remains unaltered. Thus, immediately 
behind the wave front the fluid has only a velocity normal to the wave front. In the 
region near the wave front, the normal component of the velocity, i.e. along the 
ray, is large compared to the tangential component. The wavelets, given by 


move, in the linear theory, along a ray with speed a, and the dependence of the 
amplitude on s is given by the factor [a,M6*s—M’s*]-+. This solution is now 
improved by replacing U [t—(s/a,)—7] by » and determining 7 by the condition 
that, on »=constant, (ds/dt)=a+¢,, where a is the local speed of sound and 4, is 
the velocity along the ray. 


The local speed of sound a is determined from Bernoulli’s equation 


_ 


(43) 


where y is the ratio of specific heats. If second and higher powers of ¢,/a, are 
neglected, this gives 
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In view of (41), the wavelets are determined by 


at 
[1 


and on substitution for ¢, from (39) it is seen that 


=a,~' (y+ 1) M°F (n) (46) 


After integration this becomes 


t=(s/a,)—(8a,)~* (y+ 1) M°F(n) +constant. . (47) 


[a,MB?s — M’s? 
In the linear theory the wavelets are 
U (t—s/a,—7)=constant. . : . (48) 


The constant in (47) is determined from the condition that at the body, i.e. for 
s=0, the wavelets in the linear and improved theories coincide. This gives 


U (t-—s/a,—7)+KF(y) b(s)=y, ; . (49) 


Here, and in the following formulae, the dependence of functions b, K and G 
(introduced in equation (54) ) on z-is not explicitly shown. 


It should be noticed that near the wave front, where U (t—s/a,—7) is small, the 
correction term is important. For positive acceleration b(s) will be an inverse sine 
function and for negative acceleration it will be an inverse hyperbolic sine function. 
Thus, on integration, when M’ > 0 it is seen that 


\3 
b(sy=2 (SF) sin-» (52) 
and when M’ <0, 
a,MB*\t (-—M’s\} 
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When M’ tends to zero, these expressions reduce to the corresponding ones in the 
case of steady flow. Equation (52) is not valid near s=Ma,8?/M’, for when s 
assumes this value a cusp of the wave front is reached and the solution (38) ceases 
to be valid. 


The wavelets in the undisturbed filcw propagate with speed a,, while those 
in the disturbed flow now propagate with speed a+¢,. Therefore, the latter 
(assuming that they carry an increase in pressure so that a+¢,>a,) would 
overtake the former, and this would give rise to ambiguity in the values of flow 
quantities in the region where they overlap. To prevent this, shocks must be 
introduced to separate the regions of undisturbed and disturbed flow and to dissipate 
the excess energy. As already stated, the shock speed, to the first approximation, is 
the mean of the speed of the wavelets on either side of the shock. This is easily 
seen to be true from the Rankine-Hugoniot shock conditions and is sufficient to 
determine the position of the shock. The flow quantities on the two sides of the 
shock are then given by the values on the wavelets, which meet at the shock. The 
other conditions are that the tangential component of the velocity is continuous 
and that ¢,= —a,@,. These, as noted earlier, are already satisfied in the linear theory. 


Let the position of the shock at time ¢ on a ray, characterised by the para- 
meter 7, be given by 


U (t-—s/a,—7)=-G(s). . ‘ . (54) 


Since the shock speed is the mean of the speeds of the wavelets on either side, from 
(48), (49) and (54), 


G’ (s)\=4KF (n) b’(s). (55) 


Elimination of (t—s/a,—7) from (49) and (54) gives, at the point where the wavelet 
overtakes the shock, the relation 


G (s)=KF (n) b(s)—». (56) 


Equations (55) and (56) determine G with 7 as a parameter. It is convenient to find, 
as a function of », the distance s along the ray, at which the shock cuts off the 
wavelet »=constant. If (56) is differentiated with respect to » and G’ is eliminated, 
using (55), we have 


[F? (n) b (s)] =2K—'F . . 67) 


This becomes, on integration, 


F? (n) b F (A) da, (58) 


provided that s=0 when »=0, which is true if the shock is first formed at the 
nose of the body. 
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The pressure increase at a point is given by 
P— Po= — (n) . . (59) 


which is obtained from (39), replacing £, by ». The strength of the shock, defined as 
the ratio of the pressure increase to the undisturbed pressure p, is given by 


From (58), 1 is known as a function of s, and hence (60) gives the shock strength for 
any value of s with 7 as a parameter. 


5. Details of the Flow Pattern 


In this section are collected results deduced from the theory detailed in 
Section 4. For an attached shock, formulae are obtained for the values, at the 
nose of the body, of the strength, angle of inclination with the flight path and the 
curvature. Numerical computations based on relation (60) are then given to show 
how the shock strength varies with distance and acceleration. Finally, a formula 
is deduced for the decay at large distance of a shock detached from a 
decelerating body. 


First the strength of the shock at the nose, at which s=0, will be determined to 
check that the theory reproduces the results of the steady flow problem. For small 
n, from (51), 


b(s) ~ 2s! J . (61) 


and from (40), if « is the semi-vertical angle of the nose of the body, 


Equations (58), (61) and (62) give - 


Substituting for K from (50) we get, from (62), (63) and (60), 


3 2 
Po 

for the strength of the shock at the nose. This agrees with the expression in the 
steady flow problem, since the strength of the shock at the nose depends only on 
the instantaneous velocity. 
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The angle of the shock at the nose is next determined from the equation of 
the shock (54). If the shock at the nose is inclined to the flight path at an 
angle 6, then 


From (55) Me (66) 


Hence, if lim(MG’) « 8°, 6 can be expressed approximately as 


$0 


This expression also agrees with the corresponding one for uniform velocity. It 
may be remarked that, in view of (65), since the shock cannot be inclined at an 


obtuse angle, it must be detached from the body before 6’=lim(MG’). | 


This condition can be written, after simplification, as 


m=1+(3) 


The curvature of the shock at the nose is now obtained, using the equation of 
the shock. The relations between the co-ordinates (x, r) and (7, s) of a point are 


= —f(7)—(s/M) (69) 
r=(sB/M). 
The curvature is k= =| air . (70) 
hence, in terms of s and 7, 
aM 5 +Q-B) aM & 
x=(8M)-? (71) 
4 } 
From the equation of the shock it is found that 
lim =lim 
dr M- G’ 72) 
d*s 2M’a,K*« ( 


lim =~ a9-12M4)+0 (4 | + 


129-6K‘a,22" | 


where «, is the curvature of the body at the nose. 
40 
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JSHAPE OF THE BODY 


Fig. 6. 
F-curve for body r=6{1—(1—x)*}; graphs of 5-*F(y) and 6-'r are shown. 


Equations (71) and (72) give for the curvature « of the shock 


(41 —21M*)+ O(<*) | + 


“= Bema, 
It is to be noticed from (73) that the curvature of the shock at the nose is influenced 
much more by the acceleration than by the curvature of the body at the nose. 
Cabannes obtained for the curvature of the shock at the nose an expression in the 
form 


Kee (73) 


k=g, (2, M)M’+2. (2, . (74) 


which is seen to agree in form with (73). 


Using relation (60), the shock strength is calculated for a body whose shape is 


=6, x>1 metre, 
where 8 is 0-06, so that the nose angle of the body is about 10°. This shape is one 


of those considered by Whitham"). ‘The F-curve for this body, along with its shape. 
is shown in Fig. 6 for general 5. The calculations are made for two 


_ values of Mach number, 1-5 and 1-05, and for different accelerations. For a Mach 
_ number of 1-5, the values of acceleration are 0, 6:8 and 27-2 m./sec.*. The speed 
of sound is taken to be 340 m./sec. The shock strength, plotted against the distance 


from the nose along a ray, is shown in Fig. 7. For a Mach number of 1:05, the 
shock strength is computed for accelerations 0, +3:4 and +6°8 m./sec.*. The 
results are shown in Fig. 8. It is seen that the effect of the acceleration is to 
increase the strength while that of retardation is to decrease it. The modification 
of the shock strength due to acceleration is appreciable at large distances from the 
nose and at Mach numbers near unity. In Fig. 8 the rise of the shock strength with 


_ distance in curve (5) indicates the failure of the theory near the cusp. In this case 
| the cusp is at s=1,830 m. 
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Me 1-5 
NOSE ANGLE 10°19" 
8 


P=Po 104 


4 
12 
| | 
3 H 
340 510 680 850 1020 190 ~—- [360 1530 
DISTANCE : METRES 
Fig. 7. 
(1) Uniform velocity. (2) Acceleration=6°8 metres/sec.? 


(3) Acceleration=27'2 metres/sec.? 


15 \ 
M= 1-05 
\ NOSE ANGLE 10°19' 


SS 
5 6) 
(2 
3 
340 510 680 850 1020 190 1360 1530 
DISTANCE : METRES 
Fig. 8. 


(1) Acceleration = —6°8 metres/sec.2 (2) Acceleration= —3-4 metres/sec.* 
(3) Uniform velocity. 
(4) Acceleration = +3-4 metres/sec.2 (5) Acceleration= + 6°8 metres /sec.” 
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TABLE I 
SHOCK STRENGTH IN LB./FT.? FOR A BODY OF LENGTH 40 FT., VOLUME 500 FT.®, 
NOSE ANGLE ABOUT 10° AND WITH THE SHAPE GIVEN BY (75) 


Distance from 


Mach Acceleration ‘ Steady flow Present 
number (ft./sec.?) me path Warren (Whitham) theory 
1-04 7,901 2-6 18 1-0 


1:02 —145 5,220 2:2 2°3 0-9 


For comparison, this theory is applied to determine the shock strengths in 
cases for which Warren obtained results. The values are shown in Table I. It 
may be noted that the values given by this theory are about 40 per cent. of those 
of Warren. The values given by steady flow theory are also included in Table I for 
comparison; the value increases because the distance in the second case is smaller 
than in the first case. But in the present theory, the retardation counteracts this 
effect and brings down the pressure rise. 


As the body decelerates the shock becomes detached from the nose and moves 
ahead of the body. The details of detachment are not given accurately by this theory. 
At a Mach number greater than M=1, corresponding to a ray inclined at a small 
angle to the axis, the shock becomes detached. But the ray for points on the axis is 
the one leaving the body at M=1. Thus there is a transition region in which the 
approximate theory does not apply; non-linearity effects on the position of rays 
become important here. But it is expected that away from this region the present 
theory gives satisfactory estimates of the decay of the shock. The decay of its 
strength at large distances from the body is now described. When the shock is 
detached, the lower limit of the integral in (58) cannot be zero and it has to be 
written, substituting for b (s) for M’ <0 from (53), as 


" 


where d is the initial distance of the shock, along the ray, ahead of the body and », 
specifies the wavelet on which the shock first forms. Both d and n, are dependent 
on the Mach number. But as far as the shock strength at fairly large distances is 
concerned the lower limit of the integral in (76) may be taken to be zero with 
reasonable accuracy. Further, if the nose angle is not large, d will be of the order of 
the length of the body. Hence for large s, (76) becomes approximately 


q 
F(A) da. » 


log (s/d) ~ Fe! 


(= 2M’ 


For large values of s, 7 tends to »,, the first zero of F (y), apart from »=0 (see Fig. 6). 
(If the body is of finite radius », always exists"’). 
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Therefore, from (7), Fo) ~ (= iF 


From (60) and (78), the strength of the shock for large values of s is approximately 
given by 


"0 
p Py 8a, 


[ 


) yM 0 
Py —M’/ (y+1}' s[log(s/d)}* 


(79) 


This dependence of strength on distance is typical of spherical shocks and is to 
be expected when the body drops to subsonic speeds. For the wavelets generated 
during the supersonic flight are no longer being replenished and their total energy is 
therefore limited. Thus the shock loses its similarity to a cylindrical wave spreading 
out from the axis and comes closer to a spherical wave spreading out from the region 
where the body became sonic. There are still wavelets due to the subsonic motion, 
and these overtake the shock. They will modify the foregoing result but again. 
although they have some angular variation, their amplitudes fall off with distances 
like s~', as with spherical waves. 
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Early Thoughts on the Jet Flap 


B. S. STRATFORD 


(National Gas Turbine Establishment) 


Summary: An historical sketch shows how the basic ideas of the jet flap 
were evolved. The exploratory work included the prediction that lift would 
be induced on an aerofoil by a two-dimensional jet deflected from the trailing 
edge, and that the magnitude of the total lift would be of the order of four 
times the direct jet lift. The prediction was confirmed experimentally and a 
first quantitative theory evolved. Meanwhile the hypothesis was put forward 
that, ideally, the forward thrust on the aerofoil would be independent of the 
deflection angle of the jet. 


This paper is the first of a series of papers that will describe the 
development and present state of knowledge of the jet flap. 


1. Introduction 


The present paper is the first of a series on ‘‘ The Jet Flap” to which a general 
introduction has been given in a recent Royal Aeronautical Society Lecture"’. The 
paper sketches the stages in the proving of the basic principle and in the development 
of a first quantitative theory. Later papers will be concerned with the problems 
raised by the mixing of the jet with the main stream, with certain special techniques 
that have been found necessary in experiment, and finally with theories for the more 
accurate prediction of lift. 


The jet flap principle was probably first discovered by Hagedorn and Ruden”? 
at Hanover in 1938. They observed induced lift experimentally when a two- 
dimensional jet was deflected from the trailing edge of an aerofoil, the explanation 
given being substantially that which would be put forward today. However, their 
findings were not widely published and the principle was re-discovered in 1952 by 
an entirely different approach. 


On the second occasion work had been proceeding for some time on direct jet 
lift, either for jet deflection or, with the thrust greatly boosted, for vertical take-off, 
when Constant'*’ questioned whether a gain might be achieved by using the ordinary 
but still highly powerful propulsive jet (having a thrust of the order of a third of the 
aircraft weight) to induce an external flow that would lift the aircraft. The problem 
was therefore one of magnifying the direct jet lift by a factor of say three or four. 
The first suggestion was for boundary layer control by the injection of some of the 
jet stream over a mechanical flap, this introducing the idea of the jet being arranged 
two-dimensionally along the trailing edge of the wing. The prospects looked 
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interesting, but they soon became exciting when the mechanical flap was “* removed ”’ 
and the whole propulsive jet was brought into play. 


Even after this realisation there were several months of wrangling before an 
adequate physical understanding of the system was evolved. Of the initial 
conclusions, expressed formally by the hypotheses“ quoted in Section 1.1, those 
relating to lift were developed at an early stage to yield a rough estimate of practical 
quantities. On the other hand the thrust was considered only for a much 
idealised flow. 


It should be emphasised that the present paper describes only the methods 
and ideas that were used for obtaining a first understanding; these methods would 
not necessarily be used in analysis today. 


1.1. THE HYPOTHESES 
“ The lift hypothesis 


If in a two-dimensional system a jet issues from the trailing edge of an aerofoil, 
aerodynamic lift will be induced on the aerofoil in addition to the direct vertical 
component of the jet thrust. The usual stall limitations on lift will not apply, as 
there will be a peak suction pressure at the trailing edge of the aerofoil. - Further. 
the centre of the lifting force thus created will be near the 50 per cent. chord line 
and not near the jet.” 


‘The thrust hypothesis 


Provided that the main stream speed is sufficient to prevent separation near 
the leading edge, the total forward thrust on the aerofoil will be almost independent 
of the deflection of the jet.” 


NOTATION 
main stream velocity at infinity 
main stream density 
angle of deflection of the jet at the trailing edge 


chord of jet flap aerofoil (=chord of analogous aerofoil from the 
leading edge to the beginning of the flap) 


chord of analogous aerofoil inclusive of flap 
deflection angle of analogous flap 
parameter defining size of analogous flap (as shown in Fig. 9) 


size of analogous flap relative to c 


total lift at zero incidence, i.e. jet induced lift plus direct jet 
reaction lift 
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C,, total lift coefficient based on c, at zero incidence 

Ci» pressure or induced lift coefficient based on c, at zero incidence 
| J momentum flux in jet,=reaction from jet, per unit span 

C; = J/(GpU,"c) 


G_ gain=L,/(J sin 6)=ratio of total lift at zero incidence to direct jet 
reaction lift 


d, distance of the centre of total lift aft of the mid-chord point at zero 
incidence 


a incidence of the body of the aerofoil. In most relations z has been 
put equal to zero. 


2. Induced Lift 


Suppose that the propulsive jet of an aircraft issues from a thin nozzle all 
along the trailing edge of the wing. At take-off the aircraft accelerates to some 
reasonable forward speed; then, to rise from the ground, the propulsive jet is 
deflected downwards. 


The first action of the jet on the air flowing over the wing is to mix with and 
absorb the boundary layer from the wing surface. Because of this “automatic 
boundary layer control”? an open wake or dead air region cannot exist and the 
air flow above the wing must be drawn down behind the deflected jet, as in Fig. 1(a). 
This drawing down of the upper flow creates a suction above the wing, while the 
corresponding baulking of the lower flow produces an increased pressure. Hence 
the deflection of the jet has induced a pressure lift. 


Because of the establishment of both a pressure region below the wing and a 
suction above it the approaching main stream is diverted upwards, as in Fig. 1(5), 


(3) cow 


(2) 


(1) Space to be occupied (4) High and low pressure regions established 
(2) Hence fluid drawn down (5) Fluid diverted upwards while still upstream 
(3) Therefore exhausting to “ low ” pressure 


Fig. 1. 
How the lift is induced. 
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while still upstream of the leading edge. It thus acquires an upward vertical momen- 
tum. When it turns to flow horizontally over the wing this vertical momentum causes 
an additional lift force acting on the forward surfaces. The original deflection of 
the jet has thus induced lift over both the forward and the rearward surfaces, the 
symmetry causing the centre of lift to be near the mid-chord point. 


A property of the system is that, because of the high suction pressure at the 
trailing edge corner and the automatic boundary layer control at this position, 
the boundary layer will nowhere experience a pressure rise as serious as that on 
a conventional wing. Thus the usual stall limitation will not apply. 


3. The “Times Four” Argument 


An extension of the reasoning of Section 2 shows that the magnitude of the 
lift induced by the jet is sufficient for the system to be of possible practical value. 


For simplicity only two-dimensional flow will be considered. In such a régime 
the jet, after leaving the aerofoil, must be turned until its inclination becomes zero 
at infinity downstream, for otherwise the entire field would eventually deviate from 
the horizontal and thus possess an infinite vertical momentum. The force which 
is required to turn the jet horizontal must be provided by a pressure difference 
across the jet, i.e. by an aerodynamic lift acting directly on the jet from the main 
stream. By straightforward momentum arguments this lift is such as to reduce the 
initial jet vertical momentum flux to zero, and its magnitude is therefore equal to 
the vertical component of the reaction at the nozzle (as the provision of the momen- 
tum flux is the immediate cause of the nozzle reaction). Hence, firstly, there is an 
aerodynamic lift from the main stream on the jet and, secondly, 

The aerodynamic lift on the jet is equal to the direct lift } (1) 
component of the jet reaction at the nozzle : 


The lift on the wing can now be related to the lift which has been shown to 
act on the jet. 


Two rough symmetries exist in the flow, one about an axis XX through the 
trailing edge corner and one about an axis YY through the mid-chord point (see 
Fig. 2). Due to the first symmetry the aerodynamic lift on the rear part of the 
wing, i.e. on BC, is (very roughly) equal to the aerodynamic lift on the jet, CD. 
Due to the second symmetry the lift induced on the forward part of the wing, i.e. 
on AB, is approximately equal to that on the whole, BCD, of the rear of the system; 
it is therefore equal to roughly double the aerodynamic lift induced on the jet, CD. 
Thus the total pressure lift on the wing is about three times that on the jet, or, 
from equation (1), to three times the lift from the jet reaction. 


Since the wing experiences both the pressure and also the direct jet lift, the 
ratio of the total to the direct lift is of the order of four. This ratio is referred to as 
the ‘ magnification,” or “‘ gain,” of the system. 
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Fig. 2. 
Symmetries for the “ times four ” argument. 


4, The Thrust Hypothesis 


The thrust hypothesis applies to a two-dimensional inviscid flow in which 
no mixing occurs between the main stream and the jet. A rigorous proof is 
probably difficult and will not be attempted here, but four simple arguments are 
presented. It is supposed that the jet is of high velocity so that it is sufficiently 
thin for its thickness and volume flow to be neglected; it thus becomes a laminar 
sheet having at any section a flux of momentum of magnitude, say, J (pounds per 
unit span), this sheet being inclined at the aerofoil at an angle @ to the chord line 
and being curved by the action of the main flow to become horizontal at infinity. 


(i) Consider a large “‘momentum box” around the aerofoil as in Fig. 3(a). 
When the box is sufficiently large the horizontal pressure forces will cancel and the 
net momentum flux is just the horizontal flux J (pounds) of the jet. Since the 
aerofoil provides the only solid boundary within the box the resultant pressure forces 
at this boundary must balance the horizontal flux J. Thus the aerofoil must 
experience a forward thrust of J (pounds). 


(ii) Consider thin rigid membranes AB and CD, as in Fig. 3(b), to be placed 
without disturbing the flow, on the upper and lower sides of the jet, and suppose 
that these membranes extend several chords downstream to a position where the 
jet has become almost horizontal. Let the rigid membranes be made integral with 
the aerofoil, thus forming a single body EABDC in inviscid flow with a sheet jet 
issuing sensibly horizontally from the “* second nozzle”’ BD at the second trailing 
edge. Externally, the jet will have no significant effect, as its direction is sensibly 
the same as that of both the undisturbed main stream and the rear. surfaces of the 
body; hence there will be conventional potential flow with d’Alembert’s principle 
of zero external drag applying. However, internally the body must experience the 
thrust reaction J (in order to provide the horizontal jet flow at the second nozzle 
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Fig. 3. 
The thrust hypothesis, 


BD). Thus the total thrust force on the body is J. Since by definition the mem- 
branes do not affect the flow the forward thrust J must act on the aerofoil. 


(iii) The rearward momentum flux of the jet at the nozzle is J cos 6, while 
at infinity the rearward flux is J. Thus the external flow must exert a drag on the 
jet of J(1—cos4@), as in Fig. 3(c). Using an argument similar to that of (ii), the 
total drag on the aerofoil plus jet is zero. Therefore a thrust of J (1 —cos 6) must act 
on the aerofoil. Since the direct jet thrust is J cos 6 the total thrust on the aerofoil is 


J (1 —cos #)+Jcosé=J. 


(iv) Because of the induced up-flow ahead of the aerofoil the forward 
stagnation point will be on the lower surface well aft of the leading edge. Conse- 
quently there will be a strong curvature in the flow around the leading edge, as 
sketched in Fig. 3(d), and hence a high suction pressure tending to draw the 
aerofoil forward. The suction pressures on the trailing surfaces will oppose this 
but the net effect will be to produce an induced forward thrust. The magnitude of 
the induced thrust must be J (1 — cos 4) from, say, (iii). 


5. An Experiment 


To test the hypothesis of induced lift, the pressure distribution was measured 
on the model shown in Fig. 4. 
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Fig. 4. 
The first model. 
A. One of the jet air entries. D. Jet flap. 
B. Static connections. E. Three of the static tappings 
C. Stub flap. (faintly visible within the rings). 


The model was nominally two-dimensional, so that with its span of 5 in. and 
its chord of 5} in. it reached to both side walls of the wind tunnel working section, 
which had a depth of 15 in. The width of the jet nozzle at the trailing edge of this 
unashamedly crude model varied along the span by some +40 per cent., the mean 
value being 0-0095 in. The speed of the wind tunnel could be adjusted up to 
100 ft./sec. and that of the jet, when using cold air, up to about 1,100 ft./sec., 
i.e. up to a Mach number of 1:0. The model was set at its zero lift incidence for no 
jet, this being found to be —7° when the centre line of the symmetrical “* body” of 
the aerofoil was adopted as datum. Despite the presence of the stub flap, there was 
little loss of lift from boundary layer separation just before the trailing edge, so 
that when the jet came into action substantially all the lift would be due to a direct 
influence on the main stream rather than to an effect obtained through boundary 
layer control. The deflection angle of the jet was 554° measured to the model datum 
and therefore 484° relative to the tunnel main stream. 


Figures 5 and 6 show typical pressure distributions and these suggest a sub- 
stantial pressure lift. In general the magnitude was several times greater than the 
direct lift from the jet, so that the resultant centre of lift was near the mid-chord 
position. Moreover, at a pressure lift coefficient of 2-58 (the maximum shown in 
Fig. 6) there was no sign of a stall and the only adverse pressure gradients are seen 
to be quite mild. Thus the lift hypothesis was confirmed sufficiently to justify the 
planning of a more ambitious investigation. 
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Fig. 5. 


Effect of varying tunnel speed. 
(Jet dynamic head ~ 10°12 in. Hg) 
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Fig. 6. 
Effect of varying jet speed. 
(Tunnel speed = 59 ft./sec.) 
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Fig. 7. 


Experimental lift results. 
(Lift magnification factor or lift gain). 


Figure 7 shows that the experiment confirmed also the crude “times four” 
argument, the “ gain ” being of the order of four. However, there was a significant 
variation in the gain as the strength of the jet varied. In Fig. 7 the gain is plotted 
against the “jet coefficient’? C,, =J/(4pU,°c), which is the non-dimensional 
coefficient derived from the gross thrust of the jet in the same manner as lift 
coefficient is derived from lift. When C, increased from 0-25 to 4:5 the gain 
decreased, from 6-8 to 2:3, the value four being obtained with C, close to unity. 
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This experiment had not been intended to check the thrust hypothesis (which 
was put forward after the experiment had been designed) and an attempt to show 
the existence of induced thrust met with no success. With maximum values for 
both the tunnel speed and the jet velocity, the total thrust force on the model was 
considerably less than the horizontal component of the jet reaction, which result 
could readily be ascribed to the non-uniformity of the jet nozzle width along the 
span. Verification of that hypothesis had therefore to await a later experiment. 


A final requirement was to confirm that it was the thrust coefficient rather than, 
say, the volume flow coefficient or velocity ratio which determined the lift on the 
aerofoil. Consequently pressure distributions were obtained with a mixture of 
hydrogen and nitrogen in the jet such that the density was a third of that of cold 
air and therefore approximately equal to that of a typical aircraft jet at 600°C.; for 
a given thrust this increased the velocity by a factor of 3. No difference was 
detected between the pressure distributions for the hydrogen-nitrogen jet and those, 
at the same C, values, for the cold air jet. The coefficient C, was thus established 
as the parameter determining lift when the jet angle @ is fixed. 


6. The Flap Analogy—The Name “Jet Flap” 
6.1. THE EQurvALENT AEROFOIL 


The experiment had confirmed that the ratio of the total aerofoil lift to the 
direct jet lift was of the order of four. However, although “of the order of 
four,” the ratio varied significantly with the parameter C). 


The picture that had been suggested for the flow was of the main stream air 
being drawn down behind the jet so that a region of suction pressure and therefore 
of lift became established on the rear upper surface of the aerofoil. In addition to 
this lift on the rearward surfaces, lift occurred on the forward surfaces due to the 
diversion upwards of the oncoming flow. 


Fig. 8. 
Action of the jet flap is similar to that of a mechanical flap. 


Now a similar picture holds for the flow over an aerofoil with a mechanical flap. 
(Whether the flow were potential. or real with boundary layer control, it would be 
drawn down behind the flap as in Fig. 8 and the other mechanisms just suggested 
would follow). Consequently, the downwardly extending jet sheet could be said to 
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act broadly in a manner similar to that of a mechanical flap. Although the jet 
sheet or “‘ jet flap,” as it could now be termed, reached downstream to infinity, only 
its initial portion, before it had been curved around to become almost horizontal, 
could significantly affect the flow; it would therefore induce only the same lift as 
would some specific mechanical flap of finite extent. If the strength of the jet 
increased, relative to the dynamic head of the main stream, i.e. if the value of 
C,=J/(4pU,*c) increased, the jet would penetrate further into the main stream 
and it would require a larger mechanical flap to induce the same lift as was 
induced by the jet. 


Since the lift induced by a mechanical flap in potential flow was calculable on 
standard theory, the similarity between the “‘ jet flap ’’ and the mechanical flap could 
be exploited to provide a means of calculating the pressure lift induced by the jet. 
All that was required was a knowledge of the configuration of the mechanical system, 
i.e. the angle of deflection and the chordal extent of the mechanical flap such that 
the induced lifts for the two systems would be equal. Criteria for determining this 
configuration were deduced as follows. 


In both flows the lift results mainly from the suction region above the aerofoil: 
the suction region has been pictured as being due to the drawing down of the flow 
in order to occupy and pass through the space just above the jet. The force 
required thus to draw down the flow is provided by the equal and opposite reaction 
of the flap to the lift which the flow exerts on the flap. Thus the lifts induced 
on the aerofoils would be equal if the lifts on the flaps were identical in magnitude 
and distribution. 


The requirement that the magnitudes should be equal was accepted as the first 
condition of equivalence. The second requirement could not be exactly met as it 
would require the flap shapes to be the same. The nearest that could be obtained 
to shape similarity was to make the angles of the flaps approximately equal; a more 
precise relation was considered later. The criteria of equivalence were therefore that 


(1) the aerodynamic lifts on the flaps were to be exactly equal 
(2) the flap angles were to be about equal. 


Since, from equation (1) of Section 3, the aerodynamic lift on the jet was equal 
to the direct lift component of the jet reaction, condition (1) could be simplified to 


(la) the lift on the mechanical flap must be equal to the direct lift component 
of the jet reaction, i.e. to J sin 6. 


Initially, conditions (la) and (2) were criteria for only the induced lifts being 
equal. However, the two components mentioned in condition (la) represented the 
difference between the total lifts experienced by the aerofoils and the induced or 
pressure lifts on their “‘ bodies.’’ Thus, with conditions (la) and (2) the total lifts 
also would be equal. 
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Fig. 9. 
Geometry of the analogous 
aerofoil. 
I—cosy Cy _ 2 
n= = 
l+cosf c 1+cos 


6.2. SUBSTITUTING GLAUERT’S RESULTS 

Information concerning the aerofoil with mechanical flap was obtained from 
the linearised thin aerofoil theory of Glauert‘°), the incidence of the aerofoil body 
being put equal to zero. According to Glauert, 


total lift on  . . (2) 


where c, is the chord of the analogous aerofoil inclusive of flap, 7 is the deflection 
angle of the flap, and v is a parameter defining the flap size, as shown in Fig. 9; 
while, from Fig 9, 


Since the total lifts on the aerofoils are equal, the total lift coefficient for the 
jet flap aerofoil, based on its chord c, is given from equations (2) and (4) as 


sin 2 
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The first criterion of equivalence, when used with equation (1), requires that the 
lift on the mechanical flap should be equal to the jet reaction lift, i.e. to Jsiné. 
Thus from equations (3) and (4), 


2 


Ci= sind ©) 


Equations (5) and (6) are parametric equations, relating the lift coefficient of the 
jet flap aerofoil to the jet coefficient C, through the parameter v. A direct relation 
is obtained, for at least small C,, by eliminating Y: — 


Ci =2 (27C Sin 6)! (142 + - g ). : (7a) 


which, since 2,/(27)=5-01, may for convenience be expressed as 


The absolute value of the lift may be written as 


L, =2U, (=pJcn sin 6)! (1+ a= ) 


showing that for constant jet momentum flux the lift is directly proportional to the 
main stream speed U,, except at very low speeds, and not to the square of the 
main stream speed. 


The magnification or gain, i.e. the ratio of the total lift to the direct lift, is 
obtained by dividing equation (7) by C, sin4; thus, using equation (74), 


_ sin ) 9 
(i+ 48 © ‘ (9) 
In the original parametric form the gain is 
Other relations which similarly follow from the flap analogy are: — 
= (52 sin 6 +g, * 
fsmall 
Distance of aerodynamic) 

centre aft of the aerofoil ACS 

quarter-chord point 
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d, C, sin sin \ 
and 1 120%» (13) 


where d, is the distance of the centre of total lift aft of the aerofoil mid-chord 
point at zero incidence, 


| =, sin@ ) 


where mc=the chord of the analogous flap, as in Fig. 9. 


The foregoing are the results that were obtained initially by analogy with an 
aerofoil with mechanical flap. At that time it was the gain which was primarily of 
interest, and the relationships for G (equation (9) or equations (10) and (6)) were 
therefore made independent of @ by putting »/6=1-0 for an arbitrary value of 
6 of 30°. This gave »/sin@=1-05, the resulting comparison between theory and 
experiment being shown in Fig. 7. It will be seen that the trend of the variation 
of G with C, was the same for theory and experiment, and that the magnitude of 
the experimental lift was about 70 per cent. of the theoretical. Later experiments 
gave a result nearer 90 per cent. 


This correlation with experiment pointed to the analogous flap as a useful 
device for predicting the behaviour of practical jet flap systems. The concept of 
the jet acting as a flap also gave the system its name. 
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The General Theory of “Hysteretic Damping” 


R. E. D. BISHOP 


(University Engineering Laboratory, Cambridge) 


Summary: This paper is the second of a series of three. The first“ dealt 
with the steady forced harmonic motion of a simple oscillator with one 
degree of freedom when “ hysteretic damping ” is present; by this it is meant 
that a damping force exists which is in phase with velocity but whose magni- 
tude is proportional to displacement. The present paper is devoted to the 
formulation of a general theory of small hysteretically damped vibration in 
which the existence is postulated of damping forces which act between pairs 
of points such that they are in phase with the relative velocities and their 
magnitudes are proportional to relative displacements. The theory is 
presented in general terms and use is made of Lagrange’s equations with a 
new type of “ dissipation function ”; this is comparable with that of Rayleigh 
for viscous damping. It is shown that this theory generally leads to simpler 
algebra than does that of viscous damping. Moreover, it will be shown in a 
later paper that the nature of forced and damped oscillation may be thought 
of in terms of modes more simply than is the case with the more familiar 
viscous damping. 


1. Introduction 


The purpose was discussed in an earlier paper’) of mathematical theories of 
damping in vibration analysis. It was concluded that the concept of “ hysteretic 
damping ”’ is a useful one and the theory was given for a system having one degree 
of freedom. 


A more complete mathematical theory of hysteretic damping may be based 
upon the supposition that, when a system executes steady harmonic motion, there 
acts between two points a damping force whose magnitude is proportional to their 
relative displacement and which is in phase with their relative velocity. This 
restriction of the theory to harmonic motions arises because, as shown previously. 
the definition of hysteretic damping is not wide enough to cover any but harmonic 
motions. Thus there can be no logical theory of free vibration with hysteretic 
damping until some convention is adopted for dealing with it. 


The purpose of this paper is to present, in a fairly general form, the mathe- 
matical theory of hysteretic damping for systems having more than a single degree 
of freedom. The treatment is given in terms of generalised forces and displace- 
ments and, in several respects, it resembles the classical theory of small oscillations 
with (viscous) damping. Furthermore, the results are obtained in the form of 
receptances’:*) for reasons which are explained in the text. 
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compound vibrating system 

inertia coefficient in expression for T 

sub-systems which, together, make up A 

coefficient of viscous damping 

coefficient in expression for D 

stability coefficient in expression for V 

Rayleigh’s dissipation function 

coefficient in expression for S 

flexural rigidity 

amplitude of harmonic force applied at q, 

coefficient of hysteretic damping 

spring stiffness 

length of beam 

mass 

generalised force due to presence of a hysteretic damper 
forces acting at Y and Z due to presence of a hysteretic damper 
generalised harmonic force 

component of generalised force due to hysteretic damping 
generalised co-ordinate 

hysteretic dissipation function 

kinetic energy 

potential energy 

work done during displacement in qg, only 

relative displacement of points Y and Z 


arbitrary points of a system 


displacements of Y and Z 


receptance of A 
receptances of B, C, .. . 


determinant given by equation (19) 


determinant formed from A by omitting the r column and s‘" row 
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= raf where m is mass per unit length of beam 


£,, see equation (18) 
?, amplitude of Q, 
amplitude of gq, 


frequency of oscillation. 


2. Adaptation of Lagrange’s Equations 


To generalise the theory given in a previous paper”) for a system with one 
degree of freedom, it is convenient to discuss a system having n degrees of freedom 
whose configuration is defined through the generalised co-ordinates g,, g., . . . qu. 
It will be supposed that a source of hysteretic damping is placed between two 
points Y and Z whose generalised displacements in the same direction from their 
positions in the equilibrium configuration are y and z respectively. Let the 
relative displacement be x =(y— z) so that, during harmonic oscillation, a generalised 
damping force of magnitude ihx is set up between Y and Z. 


If the damper were removed, the generalised forces exerted upon the rest of the 
system at Y and Z in the directions y and z would be, respectively, 


P, = —ih(y~—2): P.=—P,=ih(y-2). . (1) 


If y and z happen to be two of the generalised co-ordinates g, then the forces 
P, and P. will be the contributions of the damper to the corresponding generalised 
forces and will appear in Lagrange’s equations (see, for instance, Ref. 4). 


Alternatively, the co-ordinate x may be one of the q co-ordinates, in which case 
the generalised force 


will be the contribution of the damper to the corresponding generalised force at 1. 


In general, x (or y and z) will not be one of the q’s; if it is not, then 


Ox 


eq 


a Ox 


If now a virtual displacement is given to the system such that, of all the co-ordinates 
only q, is increased, then the work done by the damper force is 


bW,=P,dx= —ihxdx. . (4) 
The increment 6x in x is of magnitude 


62 The Aeronautical Quarter!s 


and 


The 


whe 
are 


dar 
sou 


fro 
fur 
for 
(0) 


so 


Fe 


= 
| 


HYSTERETIC DAMPING 


If this expression, together with that of equation (3), is introduced into (4) 
and the quantity Q,’=6W,./dq, is calculated, then it is found that 


Ox Ox Ox Ox Ox Ox 
dq, 


o/=-in| 


The contribution made by the damper to the r‘" generalised force is thus of the form 


QO,’ = +dmnGn] (r=1,2,...n) . 


where d,,, d,., and so on, are real constants. It will be seen that these constants 
are such that the order of the subscripts is immaterial; /.e. 


The results (7) and (8) have been derived under the assumption that a single 
damper exists. But they will clearly be unaltered in form by the addition of other 
sources of hysteretic damping to the system. 


It will next be shown that the forces Q,’ may be obtained by differentiation 
from a dissipation function which is somewhat similar to Rayleigh’s dissipation 
function in the theory of viscous damping’. Consider the function S which is 
formed by integrating the expression ihx with respect to x between the limits 
0 and x; thus 


When the expression (3) is substituted into (9), it is found that 
so that the generalised force Q,’ is given by 


as 


Se, (11) 


Q,’ 


_pee more, these results are not affected by the addition of more sources of 
‘etic damping to the system. 


Lagrange’s equations") may be written in the form 


34, + aq. * (r=1,2,...n) (12) 
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when they are to be used in the theory of small oscillations and when they 
incorporate the terms @D/dg, and 0S/0q, for viscous and hysteretic damping 
respectively. @Q, here represents the r‘" generalised force after removal of the 
components due to (i) viscous damping, (ii) forces having potential and 
(iii) hysteretic damping. 


3. Equations of Small Harmonic Oscillations and Receptances of 
a System having 2 Degrees of Freedom. 


The functions 7, D, V and S may be written in the form 


It follows, by use of equations (12), that the general equations of motion are 
of the form 


=P, eivt (14) 


provided that the applied forces are harmonic and of frequency © such that 
If a solution is sought of these simultaneous equations of motion such that 


then the amplitudes Y are found on substitution to be given by 


+§.¥,=9, 
(17) 
where Eve = = (c,, + i (wb,; +d,s). . (1 8) 
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The solutions for the ’’s can now be found by use of the determinant 


| S11 S12 
$21 “92 Ean 
| . 
| 
| 
Thus (r=1,2,...”) . (20) 


where the cross receptance 2,, is given by 


(21) 
A,, being that determinant which is obtained from A by omitting the r column 
and s" row. It will also be noted that, in view of the relation (18), 


It is of interest to compare these results with those of Duncan”) for the problem 
of viscous damping. 


An example will now be used to illustrate the theory of this section; it is 
convenient to take a problem which is very simple and which would scarcely 
warrant the use of Lagrange’s equations in practice. Fig. 1 shows a system which 
has two degrees of freedom. The anchored spring has, in parallel with it, a 
viscous damper and the unanchored spring has a hysteretic damper across it. 


/ 

q Felwt 

i " 


Fig. 1. 
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For the system shown 


2T = 


| (23) 
2V= +k, (q.-q +k.) qi? | | | 
“2S =ih(q. —4,)? =i + hq.” —2hq.q2] 
so that the determinant A may be written 
| (ky, +k, -M,o*)4+ +h) -(k.+ih) 
(k., + ih) +ih 
or A= {M,|M.o' — bho +k,k,} 


i {(M,+M.)h+ M.ob] -k,ob—k,h} | 


The receptances may now be written down in accordance with equation (21). 


Thus, if a harmonic force Fe“ acts at g, as shown in the figure, then 


When there is no damping (so that the quantities ¢,, are real) it will be noted 
that the determinants A and A,, are polynomials which contaifi even powers of 
only. This is still true if there is hysteretic damping present, but the expressions 
are then complex. If there is viscous damping present, however, these polynomials 
are not only complex but they also contain odd powers of ». This feature, which 
is illustrated in the given problem. leads to much heavier algebra when complicated 
systems are dealt with than is called for if only hysteretic damping is contemplated. 


This property of equations which contain hysteretic damping terms but not 
viscous damping terms is significant. That this is so may be seen from the fact 
that the general theory of conservative systems involves use of the energy 
coefficients a,, and c,, only; it follows that the effects of damping may be treated in 
general terms simply by supposing the c,, coefficients of conservative systems to 
have imaginary parts d,, added to them. 


The concept of a “complex spring stiffness’ was introduced in a previous 
paper”’. It was shown that the restoring force of a spring in parallel with a viscous 
damper may be Written in the form 


(k+ibw) x 
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where x, the spring extension, varies harmonically with frequency ». Equally, the 
restoring force of a spring in parallel with a hysteretic damper is 


(k +ih) x 


where x is again supposed to be harmonic. This suggests a convenient method of 
finding receptances when such spring-and-damper combinations are used, as will be 
shown by reference to the system of Fig. 1. 


Suppose that the two dampers in the system are removed; by elementary means 
it is then found that 


k 


2 * 


Now if k, is replaced by the complex stiffness (k,+ibw) and k, is replaced 
by (k, +ih), it is found that the previous result is obtained. 


The equations of motion (14) are given in terms of any set of generalised 
co-ordinates. Once meanings have been assigned to the co-ordinates g, a fresh set 
of co-ordinates may be taken in the form of independent linear combinations of the 
q’s in order to reduce the complexity of the equations. The “principal co-ordinates” 
of the system are a special set of this sort for which all the energy coefficients a,, 
and c,, (rs) in equation (13) vanish. However, there appears to be no reason 
why it might not sometimes be more expedient to choose a special set of co-ordinates 
in this way so as to make the energy coefficients b,, or d,, (r 4s) vanish. 


4. Some Properties of the Receptances 


Since the determinant (19) contains complex numbers, the receptances (21) 
are also complex. Now it has been pointed out by Duncan”) that complex 
receptances may be used in just the same way as real ones for the purpose of 
analysing a compound system in terms of the receptances of its sub-systems. The 
theory underlying this procedure is explained in detail in a paper by the author’? 
for real functions, and the reader will perceive that the results quoted in Ref. 3 
are immediately applicable here. 


It is convenient to illustrate this use of receptances by means of a simple 
example. Fig. 2(a) shows a cantilever which is capable of executing small, 
hysteretically damped, rotations at its root. It is stimulated by a lateral force 
Fe at its free end and it is required to find an expression for the variation of 


*This value may be taken from a table and is quoted, in fact, from entry number 8 of Table | 
of Ref, 3. 


February 1956 67 


3) 
) 
| 


R. E. D. BISHOP 


NNN 


4 HYSTERETIC 
DAMPING OF 


ROTATION 


(b) 
4s 
q (¢) 
Fig. 2. 


slope at this point of excitation. Using the co-ordinates g,, g, and q, that are shown 
in the diagram, the compound system A of Fig. 2(a) may be regarded as composed 
of the sub-systems B and C of Figs. 2(6) and 2(c); it is convenient to 
distinguish between the receptances of A, B and C by denoting them by z, 
8 and y respectively. 


Consider first the free-pinned beam which forms the sub-system B. While 
vibrating steadily with frequency , the beam will experience generalised applied 
forces at g, and g, which may be represented by ‘,,e and ®,,e. It follows 
that, if the amplitudes of g,, g, and q, are V,,, V,. and V,,, then 


=B,,%,, 


WV, =B,,%, (28) 
where £,,=/,, in virtue of equation (22). 
There is only one equation of this sort for system C, namely 
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Now in order that B and C shall fit together to form A, it is necessary to 
impose conditions of equilibrium and compatibility. First of all, equilibrium 
demands that 


®,,=F, 30 
while, for compatibility, 
Wy, =V Wyo v, 
G1) 
The imposition of these conditions upon (28) and (29) leads to the relation 
4 


Therefore the response at g, due to the excitation at gq, is given in terms of the 
receptances 8 and y of the sub-systems. 


In the problem under discussion, it is found by elementary means that 


_ cos Al sinh AJ + sin Al cosh Al 
Bis=Ba= [sin A/ cosh A/ — cos Al sinhAJ] 

_ sin Al + sinh Al 
[sinA/ cosh A/ — cosA/ sinh Al] 34 

— (cos Al + cosh Al) 
[sin AJ cosh Al —cos Al sinh AJ] 

—(1+cos Al cosh Al) 
[sin Al cosh M cos sinh Al] 


where A, / and F/ are parameters of the beam as given in the notation. If the 
damper coefficient is A, then 


(35) 


?,,, so that Y33— ih 


It follows, therefore, that z,, of equation (33) may be calculated once values have 
been assigned to the various parameters; this receptance will be complex if A is finite. 
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Results comparable with equation (33) may be tabulated for different arrange- 
ments of the sort shown in Fig. 2(c). A suitable table of this type has been 
published (see Ref. 3, Table IT). Furthermore the receptances of equations (34) 
and (35) need only be calculated once; for they too may conveniently be tabulated 
(those of equation (34) are given in Ref. 7). But it is important to notice that it may 
well be convenient to use numerical values of the sub-system receptances which have 
been found experimentally. 


5. Conclusions 


The mathematical theory has been given for systems which are subject to 
“hysteretic damping.”” The nature of the form of damping contemplated was 
discussed in a previous paper") which dealt with systems having one degree of 
freedom. 


The theory presented is somewhat similar to that of viscous damping and, in 
particular, it retains the linearity of equations of motion. It is worth noting, how- 
ever, that the algebraic expressions which appear for a system with hysteretic 
damping are simpler than those for the same system with viscous damping. The 
effect of hysteretic damping is to add imaginary parts to the stability coefficients 
Cr, Of the standard theory of conservative systems. It will be shown in a later paper 
that this leads to a simpler general theory of damped oscillation than can be had 
from use of the concept of viscous damping. 


Receptance methods are useful in the vibration analysis of damped systems and 
the general expression for any receptance is given for a system of finite freedom 
embodying both kinds of damping. The use of receptances is illustrated by means 
of an example. 
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Perturbations of Supersonic Nozzle Flows 


R. E. MEYER 


(Department of Aeronautics, University of Sydney) 


SuMMARY: Small, steady perturbations of general two-dimensional steady, 

shock-free, supersonic flows are studied and the perturbation fields of 

symmetrical nozzle flows are described in detail. The relation between errors 

in the shape of the supersonic part of the nozzle liners and the deviations 

from uniformity of the flow in the test section is given to an approximation 

sufficient for the treatment of a number of problems arising from experiment 
and in the design of nozzles. 


1. Introduction 


The study of small deviations from the design flow pattern in nozzles for 
supersonic wind tunnels is suggested by several questions arising in experimental 
work. One such problem") is the calculation, from the measured pressure distribu- 
tion in the test rhombus, of the corrections to the liner required to produce a 
uniform distribution of velocity. For this purpose it is sufficient to correct only 
a part of the liner and the problem has been solved by means of a linearised theory 
of the steady perturbations of two-dimensional simple wave flows. 


Examples of other questions arising are the effect on the velocity distribution 
in the test rhombus of a change in the stagnation pressure in order to change the 
Reynolds number, or of a slight tilting of the liners in order to change the Mach 
number slightly. For the purpose of nozzle design, means are also desirable for 
deducing, from a specified tolerance regarding uniformity of flow in the test 
thombus, the error admissible in the estimate of boundary layer corrections and, 
generally, the tolerances admissible regarding the shape of the liners. All these 
problems require an extension of the perturbation theory of Meyer and Holt"? and 
it is the purpose of this paper to present the extension to general, steady, super- 
sonic, two-dimensional, isentropic and irrotational flow of a perfect gas. 


In principle, an extension to subsonic flow is required for a treatment of 
perturbation throughout the nozzle, but this is made difficult by the change from 
elliptic to hyperbolic type of the governing equations. In fact, perturbations of 
Steady, subsonic flow are known to decay with distance from the source, especially 
when the basic flow is accelerating, and this suggests that the more. important 
perturbations are those occurring in the supersonic part of a nozzle. 
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In the following, the linearised perturbation equations of the supersonic flow 
are brought into a form (equations (13) and (14)) from which their solution can be 
obtained by a rapidly convergent process of iteration. In many cases, however, the 
full solution is not required. A rapid overall estimate of the effect is sometimes 
more desirable than a detailed prediction of its distribution and, in any case, the 
effect is often too small for its finer details to be readily measurable. The first 
approximation to the solution is then sufficient and this is presented in detail for 
nozzles with symmetrical design flow. 


The result may be summed up as follows. 


To the first approximation, errors are propagated along Mach lines of the 
design flow, so that the error at any one point of the liner is propagated along a path 
leading to the axis of symmetry, back to a point on the liner further downstream, 
thence back to the axis, and so on, until it ends at a point on the axis in the test 
rhombus. When the perturbed flow is symmetrical* the ratio of the error in Prandtl 
angle (see Section 2) at the end of this reflection path to twice the original error in 
liner slope at the beginning of the path is given by a product of reflection factors 
(equation (33)), one for each reflection at the liner. The value of A rises from 
nothing at the first reflection point for errors originating near the throat to unity at 
the end of the liner. The errors in Mach number and pressure are obtained from 
the error in Prandtl angle by conversion factors (equations (39) and (40) and Fig 3). 


In general, several points on the liner have reflection paths ending at the 
same point on the axis of the test rhombus and the total error there is obtained by 
superposition. As expected, therefore, the effect of errors in liner slope upstream 
of the simple wave region may be cancelled by that of errors farther downstream. 
When all errors in liner slope are of the same sign their effect must mount up, 
and a readily measurable effect is produced by relatively small errors. 


NOTATION 
growth factor defined by equation (4) 
h, characteristic length parameters 
m_ defined by equation (3) 
pressure 
velocity magnitude 


q 
r length of straight Mach line segment in simple wave 
s length along streamlines 


Cartesian co-ordinates in the flow plane, with x increasing along the 
axis of symmetry of the nozzle 


A __ transformation defined by equations (23) and (34a) 
A" n" repetition of transformation A 
A-"  n™ repetition of inverse of transformation A 


*The extension to asymmetrical perturbations is described in Section 6. 
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H, covariant length parameters, defined by equation (11) 
Mach number 
N _ exponent defined below equation (37) 
V defined by equation (5) 
B=0-wo 


@», Value of z in simple wave adjacent to test rhombus 


characteristic co-ordinates 


generalised characteristic co-ordinate, defined by equation (34) 


y ratio of specific heats 
€, perturbation variables, defined by equation (12) 
€,, % first approximation to €, 


6 stream direction, measured anti-clockwise from direction of x 
increasing 


6, maximum slope of supersonic liner 
A reflection factor, defined by equation (33) 
Mach angle 

o0 defined by equation (14a) 

Prandtl angle 

’ indicates small perturbation quantity 


” indicates variable of integration 


2. Perturbation Theory 


Let x, y denote Cartesian co-ordinates in the flow plane and let q, M and pu 
denote respectively the velocity magnitude, local Mach number, and Mach angle of 
the design flow. The Prandtl angle, , is defined by gdw/dq=cotyu and o=0 at 
sonic speed, and it follows from the enegy equation that 


dM |dw = M [(1+4(y- 1) M*](M?-1)°! ‘ ‘ (1) 


where y is the ratio of the specific heats. If 6 denotes the stream direction of the 
design flow, then’? 


and B=0-w . (2) 


are characteristic parameters of this flow such that z=constant on any “plus” 
Mach line, dy/dx=tan(6—,), and 8=constant on any “minus” Mach line, 
dy/dx=tan (6+ 


Let the positive direction on a Mach line be defined as that making an acute 
angle with the stream direction and let h.dz and hd denote the elements of length 
in the positive directions of the minus and plus Mach lines, respectively. 
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Then the characteristic equations of the flow can be transformed into 


0h. /0B=m(hs—h,cos 21), 


mcos 2 


and hence, if f(w)= exp | 


= {[(y—cos2)’ (sin sec p}}, . (4) 


The design flow is assumed known, that is, x, y, U and V are assumed to be known 
functions of and 


Consider now a second flow pattern which, at the point (x, y), has Mach angle, 
Prandtl angle, stream direction, and characteristic variables n+p’, w+’, 0+6, 
and 6+’, respectively, where », 0, z and are the respective values of 
these variables at the same point (x, y) in the design flow. By (2), 


a’ = 0’ + w’ and . : (8) 


Since x and y are known functions of 2 and £, these latter variables may again be 
used as independent variables. Then, if squares and products of dashed quantities 
are neglected and if conditions everywhere differ sufficiently from sonic and vacuum 
conditions for sin2 to be large compared with | 6’+,’|, the dashed variables are 
governed’) by the differential equations 


02’ /OB= - 1’), . ‘ ‘ . (10) 
where H,=(h, sin2»)~' and 


The equations (9) and (10) have been solved in Ref. 1 for the case where the 
design flow is a simple wave, that is, for regions in which either H. or H, vanish 
identically, and for the case where it is a uniform flow, so that both vanish. To 
find a general method of solution, note that, if 
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Fig. 1. 
then by equation (3) and equations (5) to (11), 

On/da= —c€, . . (14)* 
where o=m-cosec2s=c (2-8) ‘ . (14a) 


(see Fig. 1). The equations (13), (14) are similar to equations (6), (7), which have 
been treated by Meyer and Mahony in Ref. 5 and can be solved by any of the 
methods there described. 


*These equations have also been derived by Hall'4) by considering, not the difference in the 
values of (a, 8) at the same point (x, y), but the difference in position (x, y) of the two points 
where (a, 3) have the same value in the two flows. 
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For application to nozzles, however, it appears best to write (13) and (14) as 
integral equations, 


B 2 
(2, B)=E(2, B,)+ o (2” B”) aa” | ap". 
Ro 
(15) 
B 
n(z, B=n(a B)- | E02”. Bo+ | BY | da” 
Bo 
(16) 
and to solve these by iteration, starting with 


Convergence of the iteration is assured, under suitable conditions, by the results of 
Meyer and Mahony®’*; for conventional nozzles the convergence is rapid 
(see Section 4). 


The equations (13), (14) and their solution must be modified at branch 
lines’*) and in simple wave regions where the correspondence between flow plane 
and characteristic plane is not one-to-one. On a branch line H.=0, 


by (10), and on a branch line Hs=0, 


B’=constant . . . . . (20) 


by (9). In a simple wave region where H.=0 and 2=constant, 


2’= constant on any plus Mach line, (Qi) 


again by (10), since @/(hs08) represents the derivative with respect to length in the 
flow plane in the positive sense along the plus Mach lines. Moreover, by equations 
(14), (12) and (5), 


in such a simple wave region. This equation was solved by Meyer and Holt", but 
if a solution by iteration is adopted, consistent with that already proposed for 
general regions of flow, the first approximation to the solution is again relation (18). 


*A proof of convergence and bounds for the solution are also given by Hall.) 
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3. Basic Reflection in a Symmetrical Nozzle 


To apply the theory to a nozzle where both design flow and perturbation are 
symmetrical, it is sufficient to consider the upper half of the nozzle. Let 


a=A(-—f) . ‘ && 


_ be the equation in the characteristic plane* of the liner of the design flow. On the 


axis of symmetry, from which the stream direction, #, is assumed to be measured, 


a+B=0, . . (24) 


| by equation (2). Thus, if from a point B on the liner where «=<, the plus Mach line 
_ is traced in the positive sense to the axis, and thence the minus Mach line is traced in 


the positive sense to the liner (Fig. 2), then the liner is reached again at D, where 

=-—a, and 2=A(zq,). Repetition of the process leads to the point F, where 
a=A?*(z,)t. The region adjacent to, and upstream of, the test rhombus is a simple 
wave in which z=constant = z,,, say. For definiteness, assume first that z,, > A* (z,); 
then F lies upstream of the simple wave region. 


c 


Fig. 2. 


Since the perturbation equations are linear, the effect of errors in liner slope on 
the liner segment BD may be studied independently of the errors on the rest of the 
liner. Assume, therefore, that 6’ is given on the liner as a function of 2, for 
%,<a<A (z,), and that 6’=0 on the rest of the liner. The two flows compared by 
the theory are then identical for z <2, and on the Mach line BC, n=7(2,, 8)=0. 
By (18), therefore, the first approximation for » in the triangle BCD (Fig. 2) is 


and on the liner segment BD, 


a’ (a) =26’ (2), (26) 


*For better understanding of the following, reference may be made to the discussion of the 
correspondence between flow plane and characteristic plane for symmetrical nozzles in Ref. 2, 
pp. 88-92. 

+The notation A" (x) is used for the result of m successive transformations 4A (z). 
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by equations (8), (12) and (25), so that the first approximation for € in BCD is 


by equations (12) and (17), where U (z) denotes the value of U on the liner as a 
function of «. 


The first approximation for € in the triangle CDE (Fig. 2) is still given by 
relation (27). If s denotes length measured along streamlines, the streamline curva- 
ture is‘? 


d0/ds—4 sin » . (28) 

by equations (11) and (5). On the axis, 9=6’ =0, by symmetry, and hence also 


By equations (12) and (24), therefore, 


E(a, . . Gl) 


and by (18), nyo (8)=€, (— 8) & 
in CDE. 
The first approximation for » in the triangle DEF (Fig. 2) is again given by 


relation (32). By assumption, #=0 on the liner segment DF, and it follows from 
equations (8), (12), (18), (32) and (27) that 


a’ {A(z} 


where V (2) denotes the value of V on the liner as a function of z. To the first 
approximation, therefore, the errors are transmitted along Mach lines and reflected 
at axis and liner. The total effect on the error 6’ («) of transmission from liner to 
axis and back to the liner, with reflection at axis and liner, is a multiplication by 
the reflection factor, A (A (2)). 


In the absence of branch lines (so that U and V are bounded), this effect 
amounts to a reduction, without change of sign. For 2, > 0, since the theory applies 
only to the supersonic part of the flow and C lies on the sonic line when z«=0°?-°’. 
Hence, CD lies downstream of the branch line 8=0 of the nozzle and £ decreases 
in the positive direction of the plus Mach lines’, so that hs <0 and, by equations 
(4) and (5), V <0 on BC and downstream of it. On the other hand, < increases in 
the positive direction of the minus Mach lines” so that h,.>0O and U>O. It 
follows from equations (6), (7), (3) and (5), moreover, that U increases in the positive 
direction along the plus Mach lines and V decreases in the positive direction along 
minus Mach lines. By equations (33) and (29), therefore, O<A< 1. 
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It may also be noted that A—>0 as z—>0. For the Mach line CD (Fig. 2) 
approaches the branch line 6=0 and hence, 1/{V (A (z))} —>0. But the Mach 
line z=0 meets the liner downstream of the sonic line’, so U (z) remains finite. 


The process of reflection needs to be reconsidered briefly for the case 
that F (Fig. 2) lies in the simple wave region, z=constant=a,, upstream of 
the test rhombus. Assume that D lies upstream of the simple wave, i.e. 
A < % < A*(z,). To preserve continuity of notation, let a parameter be 
defined by 


= constant on any plus Mach line, 
=a,+2z/hon the axis, for z= 2,,, 


@ 
@ 


where A denotes the height of the test rhombus and z the length on the axis, 
measured downstream from the beginning of the test rhombus. Moreover, extend 
the definition of A (equation (23)) so that 


a=A(-f) . ‘ (34a) 


is the equation also of the segment in the simple wave region of the liner of the 
design flow. Since (18) is the first approximation to the solution also in the simple 
wave, and since V (a) exists also for z > zy, relation (33) then remains valid. 


Finally, consider the case that even D (Fig. 2) lies in the simple wave region. 
Equation (33) then applies only for z<,. The plus Mach line from any point on 


the liner segment in the simple wave («> 2,.) reaches the axis in the test rhombus, 
so that no reflection occurs and the first approximation to the solution is given 


directly by (21) and (26). 


It is customary to design nozzles so as to avoid discontinuities of pressure 
gradient and then’) 1/{U (z)} —>0 as z—> z,. By symmetry 1/{V (A (2))} 
and hence equations (17) and (18) must be replaced by equations (19) and (20) 
respectively. By equation (30), therefore, z’ (A (z,.))=2'(7,). The reflection factor 
A(A (z)) is thus seen to rise from nothing to unity as z rises from nothing (corre- 
sponding to a point on the liner between sonic line and throat) to 2. 


The exact values of A are most conveniently found as follows. When the 
design flow is constructed by the method of characteristics, h. and A, can be read 
off the network, by their definitions; tables of f (u), for y=1-3 and 1-4, are given by 
Meyer and Mahony”’; and A is then obtained from equations (33) and (5). When 
the design flow is constructed by the hodograph method"? 6, x and y are known as 
functions of » along the liner and U and V can be obtained from (28) and the 
analogous formula, 


dw /ds=(U-'—V~')/[4f cos 


To find V in the simple wave region, note (Ref. 1, equation (16)) that it increases 
along any of the straight Mach lines, from the upstream border of the simple wave 
to the liner, by —mr/f(u), where r is the length of the straight Mach line segment. 
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4. Accuracy of Approximation 


A detailed study of the second and higher approximations is rarely required, 
but an understanding of the type of correction they contribute to the results of the 
preceding section may be useful. From equations (15), (16), (31) and the boundary 
condition on the liner, this correction is seen to have the general form of a weighted 
mean of the first approximation. For instance, the second approximation in the 
triangle BCD (Fig. 2) is 


(2, (2), 


% 


since »(z,, 8)=0, and from the boundary condition on BC, €, is the solution of the 
integral equation 


by relations (18), (12), (27), (26) and (23). When €, (2) and €, (2) are rapidly varying 
functions, the corrections obtained from the second and higher approximations 
therefore show a tendency to even out the perturbation field predicted by the first 
approximation. In such cases the corrections may be expected to provide the major 
contribution to the full perturbations Jocally where the magnitude of the first 
approximation is particularly small, even if generally the corrections are small 
compared with the major perturbations predicted by the first approximation. On 
the other hand, if the general magnitude of the corrections is only to be compared 
with the major perturbations, then the corrections must be expected to be most 
appreciable where | « | is largest (Fig. 1) and where | U/V | reaches extreme values. 
that is, near the throat and near the test rhombus. 


The accuracy of the first approximation near the test rhombus may be further 
discussed as follows. For definiteness, take z,= 2, so that the simple wave region is 
the triangle BCD of Fig. 2 and ,=0 on its upstream border. By equation (22), (8), 
(12) and (34), and since » =constant on any straight Mach line, the distribution of 2’ 
on the liner is then the solution of the integral equation 


a’ (a) — E {2 @| adr=20 (2), . (36) 


where the integration is with respect to length on the straight Mach line 8= A-? (2), 
from the upstream border of the simple wave to the liner, the integrand being 
evaluated by the help of equation (21). Any straight Mach line of the simple wave 
may be traced, beyond the segment actually occurring in the design flow, to the 
envelope of this Mach line family and the length of the straight Mach line from liner 
to envelope is (Ref. 1, equation (16)) r,=| Ag (2) | /m. 
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Thus, if r denotes again the length of the straight Mach line segment actually 
occurring in the flow, and 2’,, the maximum of | 2’ (z)| on the liner in the simple 
wave, then by equation (36), 2’ (z) differs from its first approximation, equation (26), 
by an amount O[2’,,cr/(mr)]. In the simple wave, r=0 for z=, and (1/r))=0 
for a2=A (2), the other extreme value of z, and generally (r/r,)«1 for any con- 
ventional nozzle, since the divergence of the straight Mach lines is weak. Moreover, 
o/m remains bounded, by equations (3) and (14), as M—> oo. 


The correction contributed by the second approximation is also small near the 
throat, even if z, be chosen as small as 2°, corresponding to M = 1-1 on the liner, 
i.e. a liner point very close to the throat*. For on the one hand, the liner segment 
BD (Fig. 2) then crosses the branch line 8=0 and remains close to it™, so that 
| U (z)/V (a) | in equation (35) is small. On the other hand, although «—> 0 as 
M —> 1, the maximum value occurring is « {2° + A~'(2°)} ~5 and o falls rapidly 
with increasing argument to vanish for M=2(3-— y)~!; moreover, 


=A(a,) 


where tan 6,, is the maximum slope of the liner, which is rather small. 


Generally, therefore, two facts combine to make the difference between first 
and second approximations rather small. Firstly, this difference has the form of a 
sum of integrals, each over an interval bounded by @,,. Secondly, the integrals are 
multiplied by factors which are small near the throat, because of the branch line of 
the symmetrical design flow, and also further downstream, because of the weakening 
divergence of the Mach lines, due to the progressive decrease in magnitude of the 
pressure gradient with distance from the throat in conventional nozzle designs. 
Everywhere, therefore, this difference is rather small compared with the major 
perturbations predicted by the first approximation and, in the following, considera- 
tion will be restricted to the first approximation. 


5. Resultant Deviations 


To any point on the axis of the test rhombus there corresponds a point on the 
liner, in the simple wave, where a has the same value. By equation (21), 2’ also has 
the same value at both points. By equation (33), therefore, the first approximation 
for the total error 2’ on the axis of the test rhombus, due to all errors & on the 
supersonic part of the liner, is 


N n-1 


a! (z)=20' (@)+2 ILA{A-*(@)} {A-"(@)}, 


n—1 k=0 


where N (2) is the greatest integer such that A~‘(z)>0. If both design flow and 
perturbation are symmetrical, 


w’ (a)= 2’ (2) ‘ . (38) 
on the axis, by equations (8) and (30), and 


*None the less, sin 2u>0°7. 
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by equation (1), where M is the final Mach number of the design flow. Similarly, 
by the energy equation and isentropic relations, | M 


= — (M? — w’ (40) line 
where p is the final pressure of the design flow. The two conversion factors are 
shown in Fig. 3, for y= 1:4. The deviations at other points of the test rhombus are inct 
obtained from those on the axis as described by Meyer and Holt’. of : 


As an example, consider the nozzle of Ref. 2, where M=2:0, 2,=26° and A 
A~' (ay) = 3°, so that N=1 for most values of z in the test rhombus, but N =2 for ens 
the highest values occurring. From the figure of Ref. 2, A(z.) 0-1. Thus if 6’ (z) test 
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/ 
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were constant along the liner, the resultant deviation in Mach number would rise 
along the axis of the test rhombus from 4:56’ to 8:5, approximately, and the 
relative deviation in pressure, from 76’ to 136’, approximately. If a tolerance, 
| M’| <«, were specified for the uniformity of the Mach number distribution in the 
test rhombus, it could be ensured, for instance, by computation of the distribution of 
liner slope (including boundary layer correction) to within +¢/8°5. 


The resultant deviations caused by given errors are somewhat lower in nozzles 
for Mach numbers less than 2, but they increase progressively as the Mach number 
increases, for not only the conversion factors increase (Fig. 3), but also the values 
of z’ obtained from (37). For instance, for the nozzle of Ref. 8, M=2-41, 2, =37°, 
(aw) 15°, A-? (a) 1°, so that N=2 throughout, and A(z.) ~1/3. Thus if 
#=constant, M’ rises from 6°40 to 116’, approximately, along the axis of the 
test rhombus. 


The distribution #=constant represents approximately the effect of a change 
in stagnation pressure", if the deviations in liner slope upstream of the throat and 
the small change in throat width thereby caused are neglected*. If the nozzle of 
Ref. 8 were corrected according to the estimates there given, but run at a stagnation 
pressure five times higher than assumed there, to obtain a Reynolds number five 
times higher, the actual rate of growth of the displacement thickness would be 
below that corrected for by a factor 5-'/°, i.e. by about 8 x 10-* on the average. 
The Mach number would then increase along the axis of the test rhombus by 
approximately 0-004 and the mean Mach number would be above the design value 
by about 0-007. Usually, therefore, the effect of changes in stagnation pressure on 
the velocity distribution in the test rhombus should be too small to affect experi- 
ments, but it may be noticeable in calibration tests. 


6. Asymmetrical Perturbations 


If the design flow is symmetrical, but the errors in liner slope are not, equation 
(37) still describes the contribution of the upper liner to the deviations on the axis 
of the test rhombus, but equation (30) does not hold there, nor equation (38), for 
the lower liner contributes deviations 6’(a). These can be calculated from the 
results already given by symmetry, and «’ (z)=4 { 2’ (z)— 8’ (2)}, and the deviations 
of stream direction on the axis of the test rhombus are given by 4{2’(z)+/’ (2)}, 
by equation (8). 


A case of more importance is that where the nozzle is designed to have one 
plane liner, equivalent to the axis of symmetry. The effectiveiiner, however, is not 
quite flat, due to the boundary layer growth on it. In that case, equation (37) still 


describes the contribution of the curved upper liner, but the nearly flat, lower liner 
may contribute further deviations. If 4, (a) denotes the deviations of liner slope on 
it, then those occurring between the sonic line and the test rhombus contribute 


deviations 


*For simplicity, the boundary layers on the side walls are also neglected. 
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on the axis of the test rhombus, where N is defined as for equation (37), and the 
reflection factor A is still given by equation (33), forO<k <n-—2, but 


= —Vi (A-"(@} /V 


where V,(z) denotes the distribution of V on the plane liner as a function of «. 
This new reflection coefficient represents only the effect of transmission from axis to 
liner and reflection there; it is negative, since V <0, and tends to —1 both as 
A~"(a)—>0 and as A~"(z)—> 2, by equation (20), but otherwise is less than 
unity in magnitude, by equation (7). Deviations 6’,(@) occurring on the nearly flat 
liner in the test rhombus, on the other hand, contribute no deviations 2’ and hence, 
by equation (8), their contribution to the deviations in Prandtl angle is —6’,(a). 


Altogether, therefore, w’ (a) =a’ (2) + 2’, (2) 6, (2), 


from which the deviations in Mach number and pressure are found by equations 
(39) and (40). 


The theory can also be applied to nozzles where the design flow is not 
symmetrical, e.g. sliding-block nozzles, but in such cases the reflection paths must be 
traced from liner to liner, and back, and it is difficult to obtain general results 
concerning the reflection factors. 
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